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1 .  Backoround:  Review  of  Basic  Concepts 

This  final  report  illustrates  the  results  of  a  three-years  research 
program  on  Quantum  Chaos.  When  our  program  was  started  this  subject 
was  still  in  a  native  stage,  dominated  by  the  pioneering  work  of  a 
restricted  number  of  scientists.  By  now,  instead,  it  has  grown  into  a 
most  active  field  of  research,  and  keeps  drawing  increasing  attention  by 
scientists  working  in  widely  different  areas,  ranging  from  pure 
theoretical  to  experimental  ones.  International  conferences  and 
workshops  on  Quantum  Chaos  have  been  held,  and  others  are  being 
announced,  witnessing  the  liveliness  and  the  interest  of  the  topic. 

Curiously  enough,  despite  this  great  development  the  very  definition 
of  Quantum  Chaos  is  still  a  controversial  point.  A  very  neutral  statement 
is  that  quantum  chaos  is  concerned  with  the  properties  of  quantum 
systems  that  are  chaotic  in  the  classical  limit.  In  order  to  get  a  more 
precise  formulation  of  the  problems  involved  and  to  give  a  clear 
exposition  of  our  contributions,  we  shall  first  review  some  basic 
concepts  and  keywords. 


1.1  Classical  Chaos. 

Chaos  is  today  a  protagonist  in  all  problems  that  can  be  treated  by 
classical  mechanics.  The  appearance  of  chaos  is  connected  with  an 
extreme  type  of  instability  that  often  appears  In  nonlinear  classical 
systems.  Chaotic  systems  are  in  many  senses  the  very  opposite  of 
integrable  systems.  The  orbits  in  phase  space  ot  one  integrable  system 
are  bound  to  smooth  surfaces  whose  dimension  equals  the  number  of 
freedoms,  called  invariant  ton.  For  this  reason,  integrable  systems  are 


veru  itatile;  a  itnall  error  in  speciluinq  the  initial  data  involves  an  error 
in  predicting  the  state  at  a  later  time  that  grows  linearlg  with  time,  so 
that  finite-precision  algorithms  allow  for  a  reliable  computation  of 
orbits.  Integrable  systems  are  well  known,  being  the  sole  class  of 
nonlinear  systems  that  are  amenable  to  analytical  solution. 

Generic  Hamiltonian  systems  can  be  regarded  as  perturbations  of 
integrable  systems  -  i.e.,  their  Hamiltonian  is  obtained  by  adding  some 
correction  to  an  integrable  Hamiltonian.  It  is  now  known  that  in  this 
case  some  orbits  are  still  confined  (and  densely  till)  on  some  invariant 
surfaces  that  can  be  looked  upon  as  distorted  remnants  of  the  supply  of 
invariant  ton  that  is  available  in  the  unperturbed  case.  The  celebrated 
KAfl  theorem  characterizes  this  surviving  set  of  'regular'  orbits,  that 
becomes  narrower  as  the  perturbation  strength-or  just  the  energy,  in  the 
conservative  case-increases. 

Chaos  became  a  keyword  in  classical  mechanics  when  the  behaviour  of 
the  other  orbits  -  i.e.,  those  that  are  no  longer  confined  to  smooth 
''tori''-was  studied  in  some  detail,  indeed,  it  turned  out  that  these  orbits 
are  usually  exponentially  unstable  -  i.e..  errors  In  the  initial  data 
propagates  exponentially  in  time.  This  makes  hopeless  the  prediction  of 
individual  orbits  by  finite  precision  algorithms.  If  the  state  of  the 
system  moving  along  a  typical  chaotic  orbit  is  recorded  with  a  finite 
precision  at  a  prescribed  sequence  of  instants,  one  gets  a  string  of  data 
that  is  indistinguishable  from  a  random  string  -  i.e.,  from  a  string  of 
numbers  generated  by  a  random  device.  Analysis  of  single  orbits  is 
therefore  meaningless;  instead,  the  behaviour  of  ensembles  of  orbits 
turns  out  to  be  describable,  in  statistical  terms,  and  it  is  found  that  the 
ensemble  evolution  is,  in  a  suitable  approximation,  a  diffusion  process. 
Whereas  regular  orbits  lie  on  invariant  ton  and  are  therefore  confined 
forever  to  a  negligible  region  of  phase  space  chaotic  orbits  have  a 
tendency  to  explore  the  whole  available  phase  space  in  a  diffusive  v/ay. 
Bounds  on  this  pervasive  behaviour  are  however  posed  by  tne  survival  of 
some  tori,  that  can  substantially  slow  down  or  even  stop  the  diffusion. 


I..L  riodei  sysferris  -  Euliards  and  the  iicied  Rotator 

The  onset  of  chaos  in  Hamiltonian  systems  is  in  most  case; 
unambiguuusly  detected  by  numerical  simulation,  but  the  tiieoretn;  al 
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jfialysis  of  this  process  is  still  far  froro  being  cornplete.  For  tins 
reason,  in  order  to  understand  at  least  the  qualitative  features  of 
chaotic  dynamics  it  has  proved  very  useful  to  consider  simplified 
models.  Not  being  modelled  after  any  actual  phenomenological  situation, 
they  are  abstract  constructions;  yet,  they  exhibit  the  essential  features 
of  chaos,  bare  of  complications  that  are  usually  superimposed  in 
realistic  cases.  A  first  example  is  provided  by  billiards,  i.e.,  particles 
bouncing  elastically  inside  suitable  plane  regions.  Rigorous  mathematical 
results  have  shown  that,  if  the  boundary  of  the  billiard  is  suitably 
shaped,  then  the  simple  Hamiltonian  system  so  defined  displays  strongly 
chaotic  properties.  Unlike  generic  Hamiltonian  systems  in  which  both 
chaotic  and  regular  orbits  coexist,  the  latter  becoming  negligible  only  by 
suitably  increasing  the  energy,  billiards  are  fully  chaotic  -  i.e.,  regular 
orbits  have  zero  measure  -  at  any  energy. 

As  we  shall  emphasize  later,  the  most  Impressive  phenomenological 
manifestations  of  chaos  are  offered  by  Hamiltonian  systems  subject  to 
external  perturbations  periodic  in  time.  Within  this  class,  one  model 
system  dominates  the  scene:  the  6-kicked  rotator,  also  called  standard, 
or  Chirikov's  map,  by  more  mathematically  oriented  authors.  This  is  the 
i-dim.  Hamiltonian  system  described  by  the  time-dependent  Hamiltonian 

H=  P^/2  +  kcos  0  2  S(t  -  nT);  o  <  <  2Tf 
n 

Physically,  it  is  a  pendulum  whose  weight  is  'turned  up'  in  an  impulsive, 
6-like  way,  at  regularly  operated  instants  of  time.  It  can  obviously  be 
looked  upon,  as  a  perturbation  of  a  free  rotator  (k=0)  and  indeed,  the 
whole  scenery  of  KAM  theory  is  here  realized.  The  object  of  interest 
here  is  the  discrete-time  orbit  p(nT),  0(nT),  which,  in  the  unperturbed 
case  is  generically  bound  and  densely  fills  a  regular  curve.  As  k  is 
increased,  the  set  of  'regular'  orbits  shrinks  and,  for  kT>l.  almost  all 
orbits  become  chaotic.  A  remarkable  indication  of  the  stochattc  nature 
of  the  motion  in  this  case  is  offered  by  the  behaviour  of  the  kinetic 
energy  T(nT)=  <(P-(nT)/2>  averaged  over  the  phase  0;  this  grows  in  a 
diffusive  way  ~  On  with  D  %  (kT)^/2. 


-r 
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1.?  Pheiiufiit'i'iulijQical  Rf'levanc £■  of  Llassical  Chao:- 

The  appearance  of  irregular  orbits  may  dramatically  change  the 
phenomenological  behaviour  of  a  given  system.  Indeed,  a  major  impulse 
to  the  study  of  chaotic  dynamics  was  given  by  the  necessity  of 
predicting  the  precise  conditions  under  which  dangerous  instabilities 
develop  in  physical  situations  that  are  amenable  to  a  classical 
description,  such  as,  e.g.,  beams  of  accelerated  particles  or  confined 
plasmas.  A  vivid  illustration  of  the  effect  of  chaos  on  real  physics  is 
provided  by  the  very  problem  that  was  the  most  important  part  of  the 
research  described  here  :  an  hydrogen  atom  in  an  external  microwave 
Held.  Deferring  a  detailed  analysis  to  a  later  section,  we  shall  give 
here  a  qualitative  explanation  why  chaos  is  essential  in  the  classical 
behaviour  of  such  an  atom.  The  unperturbed  atom  is  obviously  an 
integrable  system.  As  soon  as  the  external  perturbation  is  turned  on  (the 
microwave  field),  some  of  the  orbits  become  irregular  and  start 
wandering  away.  For  not  too  large  perturbation  strength,  residual 
invariant  tori  persist  that  prevent  this  diffusion  from  leaving  the  'bound' 
state  subspace  of  the  whole  phase  space.  However,  for  any  given  initial 
(bound)  state  of  the  atom,  a  value  of  the  external  field  will  be  found,  for 
which  the  orbit  leaving  from  this  st  .e  will  no  longer  meet  any  such 
Impediment  and  will  therefore  diffuse  away  until  the  atom  ionizes.  This 
simple  picture  leads  to  predict  that  for  atoms  prepared  in  a  fixed  initial 
state  under  a  field  of  fixed  frequency  there  is  a  threshold  in  the  field 
strength,  across  which  the  ionization  rate  changes  abruptly  due  to  the 
onset  of  unbounded  diffusion. 


l.d  Quantum  Limitations  to  Classical  Chaos  in  Conservative  Systems. 
Level  Statistics. 


Why  should  a  quantum  physicist  care  about  classical  chaos'''  In  the 
light  of  the  above  sketched  picture  of  the  microwave  ionization  of 
hydrogen  atoms,  this  question  has  at  least  one  immediate  answer 
because  just  the  lailure  of  certain  classical  prediction::,  about  the 
hydrogen  .atoms  'was  one  starting  point  in  the  edification  of  quantum 
mechanics.  Nevertheless,  on  account  of  the  correspondence  principle,  one 
can  predict  that  by  going  up  to  sufficiently  higli  quantum  numbers 
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quantum  and  classical  predictions  will  agree.  It  was  however  apparent 
since  the  beginning  of  studies  on  chaotic  dgnamics  that  the 
correspondence  between  a  quantum  system  and  its  classical  limit  may 
have  some  subtle  aspect  when  the  latter  is  chaotic.  For  example,  a 
chaotic  conservative  classical  system  with  a  bounded  configuration 
space  (such  as  a  billiard)  display  a  highly  non-recurrent  behaviour:  by 
this  we  mean  that,  even  though  it  comes  infinitely  often  arbitrarily 
close  to  its  initial  state  (Poincare  recurrence),  nevertheless  there  is  no 
upper  bound  on  the  return  times,  that  will  take  arbitrarily  large  values 
as  the  system  moves  along  a  given  orbit.  The  quantum  analog  of  any  such 
system  will  have  a  pure  point  energy  spectrum;  the  quantum  evolution 
will  then  be  quasi-periodic,  hence  recurrent.  On  the  other  hand,  the 
correspondence  principle  requires  that  the  recurrent  quantum  behaviour 
goes  into  the  highly  non  recurrent  classical  one  as  h  ->  0! 

An  obvious  way  out  to  this  apparent  contradiction  is  that  recurrence 
(or  non-recurrence)  is  a  long-time  property,  and  the  classical  limit  fi->o 
can  be  non  uniform  in  time.  In  other  words,  for  a  given  quantum  state 
which  is  quasi-classical  one  can  expect  at  best  that  quantum  dynamics 
looks  like  the  classical  over  just  a  finite  time  scale  r(1i)  so  that 
r(1i)->oo  as  fi->o.  One  has  then  a  sort  of  chaos  confined  to  a  bounded 
time  scale,  so  that  "Transient  Chaos"  or  "Pseudochaos"  appear  more 
appropriate  expressions. 

Despite  this  severe  limitation  placed  by  quantum  mechanics  on 
classical  chaos  in  the  conservative  case  -  the  nonconservative  case  will 
be  discussed  below  -  there  are  still  reasons  why  quantired  conservative 
chaotic  systems  should  be  made  the  object  of  a  careful  investigation. 
Integrable  systems  played  a  central  role  in  the  early  stages  of  quantum 
mechanics.  Indeed,  the  Bohr-5ommerfeld  quantization  procedure  is  only 
applicable  to  integrable  systems. 

Instead,  the  present-day  form  of  quantum  mechanics  does  no  longer 
suffer  from  such  a  limitation,  and  we  can  safely  quantize  any 
conservative  Hamiltonian  system,  regardless  of  its  integrable  or  chaotic 
nature.  Nevertheless,  semiclassical  quantization  rules  still  play  an 
important  role,  in  that  they  are  often  the  only  reasonable  way  to  get 
quantitative  informations  about  the  higher  part  of  the  energy  spectrum. 
These  rules  are  today  rigorously  deduced  from  tlv  Schroedinger  equation 
in  the  form  of  the  EWBK  rules.  Again,  this  theory  works  only  under  the 
explicit  assumption  that  the  system  at  hand  has  an  integrable  classical 


limit!  A  rijtuiol  tiroblt'iTi  is  then  whether  semic l.assicol  methods  can  be 
used  also  m  the  non-integrable  case  in  order  to  get  infromation  about 
the  serniclassical  part  of  the  spectrum. 

(t  was  suggested  by  Berry  111  that  the  integrable  or  non-integrable 
character  of  classical  dynamics  is  mirrored  by  the  statistical  type  of 
the  corresponding  energy  spectra.  Given  a  spectrum  whatsoever,  one  can 
define  the  associated  level  spacing  distribution  function  P(s,E)  as 
follows.  Given  a  A>0,  one  takes  only  those  eigenvalues  Ep(fi)  that  lie 

between  E-A  and  E+A  and  forms  the  string  s^  of  the  spacings  Ef,+  ,  -  E^^, 
measured  in  units  ot  their  average  value.  These  numbers  s,-,  'will  be 

distributed  according  to  an  hystogram  N(s.fi.E)  giving,  for  any  s  the 
relative  frequency  of  occurence,  within  the  string,  of  spacings  between  s 
and  s+ds.  Letting  then  ti->o,  Nts,fi,E)  -would  tend  to  a  limit  P(s,E). 

On  account  of  numerical  computations  and  analytical  arguments,  it  is 
currently  assumed  that.  If  the  system  has  an  integrable  classical  limit, 
P(s,E)  is  given  by  the  Poisson  distribution  e’’’.  Instead,  if  the  system  has 
a  chaotic  classical  limit,  P(s)  should  be  the  same  as  the  Wigner 
distribution  valid  for  eigenvalues  of  random  matrices  in  the  Gaussian 
orthogonal  ensemble,  that  is  given  almost  exactly  by 

Pis)  :  (n/2)S  exp(-ns^/4) 


The  most  obvious  difference  between  these  two  forms  of  P(s)  is  that, 
in  the  latter,  one  has  level  repulsion  -  i.e.,  the  frequency  of  small 
spacings  is  vanishingly  small. 


1.5  Quantum  Chaos  for  Periodically  Perturbed  Hamiltonian  Systems. 

In  tlie  section  4  above  a  reason  was  pointed  out,  w'hy  the  long  -time 
behaviour  of  quantum  and  classical  conservative  systems  in  the  presence 
of  chaos  cannot  be  expected  to  be  the  same.  However,  the  physically 
most  interesting  cases  such  as  the  H-atorn  of  sec.  1.?.  do  not  belong  to 
this  class,  because  their  Hamiltonian  is  periodic  in  time. 

In  that  case,  quantum  dgnamifs  is  no  longer  denuod  by  an  energy 
spectrum,  but  rtther  by  the  quasi-energy  '.q.e.i  spectriim.  'while  in  the 
conservative  case  the  energy  spectrum  was  ino'wn  to  be  pure  point,  there 


li  now  iKi  reasun  to  expect  that  the  q.e.  spectrum  should  be  pure  point  m 
all  cases;  the  question  of  persistence  of  chaos  or  of  some  of  its  effects 
in  quantum  mechanics  in  this  nonconservative  case  must  then  be  started 
atresh. 

The  first  results  were  obtained  bg  Casati  et  al.  14|  who  subjected  the 
quantum  version  of  the  S  kicked  rotor  sketched  m  sec.  1.2  to  a  numerical 
analysis, with  the  following  remarkable  results.  It  will  be  recalled  that 
the  character  of  the  classical  motion  is  defined  by  the  quantity  K-kT.  In 
the  quantum  case,  the  dynamics  depends  on  two  parameters  k^=kfi-’  and 

Tq-Tfi.  If  Tq  is  a  rational  multiple  of  dfl  dn  m/n  ,  rn,  n  integers)  the 
quantum  rotor  behaves  resonantly,  and  its  average  kinetic  energy 
E(t)=<'f(t)  I  -  (fi^/2)d^/d0^  I  t(t)>  increases  with  time  according  to  an 
asymptotically  t^  law.  This  resonance  has  no  classical  analog  and  is  not 
therefore  related  to  classical  chaos.  Therefore,  the  search  for  quantum 
chaos  must  resort  to  irrational  values  of  (4n)''T,j.  The  response  of 

computer  experiments  was  that  E(t)  follows  the  classical  diffusive  law 
E(t)at  only  up  to  some  time  t^,  after  which  it  enters  a  steady 

oscillatory  state.  In  other  words,  according  to  numerical  experiments 
one  has  a  quantum  suppression  of  chaotic  diffusion. 


1.6  Relevance  of  Classical  Anallsys  to  the  Problem  of  Microwave 
Ionization 

The  other  nonconservative  system  discussed  in  this  report  -  the 
H-atorn  in  a  microwave  field  -  had  never  been  subjected  to  quantum 
anallsys,  neither  theoretical  nor  computational,  up  to  the  start  of  our 
program.  Nevertheless,  it  presented  a  unique  occasion  to  check  the 
possible  survival  of  chaotic  effects  in  quantum  mechanics,  tiecause 
laboratory  experiments  on  micro'wave  ionization  for  highly  excited 
H-atoms  had  been  performed  since  1974  by  Bayfield  and  Koch.  [12]  who 
succeeded  in  exposing  a  relevant  ionization  for  low  fields  and 
frequencies  well  below  the  l -photon  threshold.  No  satisiartory  quantum 
analysis  was  available  tor  this  phenomenon.  Alter  tne  won  ot  jensen 
119.20]  the  idea  took  place,  that  quantitati'-'e  predictions  in  agreement 
vvith  experiments  might  be  provided  by  a  classical  theory;  indeed. 
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uverjll  ijtist'dctory  response.  Accordiny  to  such  results,  quantum 
mechanics  does  indeed  loilow  classical  mechanics  so  cioseiu,  as  to 
rpproduce  even  chaotic  ettects. 


2.  statement  of  Problems 

The  rapid  overview  of  the  state  of  the  art  we  gave  in  the  previous 
section  should  have  put  m  evidence  that  current  prohlerns  in  Quantum 
Chaos  Delong  to  two  main  areas:  those  involving  the  spectral  properties 
of  conservative  systems  and  those  investigating  the  dynamics  of 
externally  perturbed  .systems. 

In  this  report  we  shall  present  a  number  of  contributions  we  have 
given  to  both  fields. 

2. 1  The  conservative  case 

Our  starting  point  here  was  Berry’s  hypothesis,  that  the  statistics  of 
energy  levels  should  be  essentially  different  in  the  integrable  and  in  the 
chaotic  case.  Numerical  work  by  Bohigas  and  other  authors  [1]  supports 
the  generally  accepted  conclusion,  that  the  level  statistics  for  a 
quantum  system  which  is  chaotic  in  the  classical  limit  is  well  described 
by  random  matrix  theory:  in  particular,  under  the  assumption  of  time 
reversal  invariance  one  can  assume  that  fluctuation  properties  in  the 
spectrum  will  obey  tne  same  statistics  as  in  the  case  ot  random 
matrices  in  the  so  called  Gaussian  Orthogonal  ensemble.  This  contention 
was  tested,  by  analyzing  "Ist  order"  statistics-  i.e.,  statistical 
properties  involving  pair  correlations  between  different  levels,  such  as 
the  spacing  distribution  Pts),  and  also  higher  cder  statistics,  involving 
multiple  correlations,  such  as  the  6-3  statistics  to  be  defined  later. 
Therefore,  chaotic  systems  display  a  remarkable  universality  in  their 
■-.pectr.al  statistics. 

I-;,  there  any  universal  statistics  tor  the  integrable  case''’  Early  resuU.s 
suggested  an  aiiirmati'-'e  ans'»ver;  indeed.  Berry  gave  an  argument 
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distribution  of  spjcmqs.  On  the  grounds  of  this  orgument  it  has  been 
conjectured  that  the  fluctuation  properties  of  the  'integrable'  spectrum 
should  be  the  same  as  tor  a  Poisson  process.  According  to  this 
conjecture,  the  energy  levels  of  an  integrable  system,  despite  their 
predictability  by  more  or  less  simple  semiclassical  rules,  should  make 
up  3  completely  random  sequence! 

Even  though  this  conjecture  was  supported  by  numerical  results 
insofar  as  the  level  spacing  distribution  was  involved,  there  was  neither 
theoretical  nor  numerical  evidence  that  also  higher  order  statistics 
should  be  Poisson-like.  Therefore,  we  posed  the  following  problem: 
fo  analyze  higher  order  statistical  properties  -  in  particular,  the  6-5 
statistics  -  of  the  sequence  of  energy  levels  of  an  integrable  system,  in 
order  to  check  whether  Poisson  statistics  is  still  obeyed. 


2.2  The  Time  Dependent  Case 

Our  current  understanding  of  the  interaction  between  matter  and 
radiation  is  provided  by  quantum  mechanics,  the  development  ot  which 
was  indeed  prompted  by  the  incapability  of  classical  mechanics  to 
account  for  the  response  of  microsystems  to  radiation  fields.  In  simple 
cases,  a  sufficient  approximation  for  the  behaviour  of  an  atom  or  a 
molecule  to  an  external  electromagnetic  field  con  be  obtained  by 
considering  the  latter  as  a  classical  field;  in  these  cases,  the  quantum 
dynamics  is  described  by  a  time  dependent  Schroedinger  equation. 

When  the  intensity  of  the  fields  is  sufficiently  small,  this  can  be 
handled  by  time-dependent  perturbation  theory,  that  leads  to  describe 
the  interaction  with  radiation  in  terms  of  multiphoton  processes. 
However,  the  present  state  of  physical  research  often  confronts  us  with 
situations  in  which  it  is  important  to  predict  the  response  of  atoms  or 
molecules  to  intense  fields;  moreover,  experimental  results  indicate  that 
the  dynamics  in  such  cases  may  be  qualitatively  very  different  than  one 
would  guess  on  the  grounds  of  perturbatively  based  intuitions.  A  typical 
example  is  Bayfield  and  loch's  experirnerd,  which  exposed  strong 
ionization  in  a  situation  in  winch  '  inc)  phorons  would  be  required. 

In  the  absence  of  a  method  whatsoever,  by  vvhich  the  Schroedinger 
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required,  the  best  one  oon  do  m  order  to  get  some  theoretical  Indication 
IS  just  coming  bact  to  classical  mechanics,  and  indeed  in  physically 
relevant  cases  such  as  microwave  ionization  oi  hydrogen  atom  in  an 
external  field  shows  that  a  chaotic  threshold  exist.  For  field  strength 
exceeding  this  threshold,  a  qualitative  change  in  the  dynamics  occurs, 
that  leads  to  intense  ionization  due  to  a  diffusive-like  motion  of  the 
electron  in  the  external  field.  Numerical  simulation  of  the  classical 
model  IS  easily  feasible,  and  a  partial  comparison  of  numerical  data 
gotten  in  this  way  with  results  of  experiments  on  real  atoms  has  shown 
a  certain  degree  or  agreement.  This  fact  seems  to  indicate  that  in  a 
semiclassical  regime,  the  predictions  ot  classical  chaotic  dynamics  are 
essentially  respected  by  quantum  dynamics,  ohall  we  conclude  that 
classical  chaos  in  time-dependent  problems  is  basically  surviving 
quantization?  Certainly  not;  and  the  reason  is  that  we  know  of  at  least 
one  quantum  system  -  the  .S-kicked  rotator  -  in  which  the  fully 
developed  chaos  of  the  classical  model  is  completely  suppressed.  What 
are  then  the  reasons  why  quantum  mechanics  reacts  so  differently  to 
classical  chaos  in  these  two  cases?  Answering  this  question  is  an 
essential  task,  because  the  existence  of  a  quantum  regime  lying  beyond 
the  perturbative  regime  in  which  the  motion  has  some  ’diffusive' 
features  would  open  a  new  field  of  immense  potential  interest.  We 
fhrefore  need  some  theory  that  may  be  able  to  set  precise  quantitative 
conditions  for  the  applicability  of  semiclassical  approximations  when 
the  classical  motion  is  chaotic.  The  following  program  should  then  be 
fulfilled: 

1)  To  build  a  simple  model  for  a  H  atom  in  a  microwave  field,  that  is 
amenable  to  both  classical  and  quantum  solution  by  computer  simulation 

2)  to  compare  the  classical  and  quantum  motion  in  different  parameter 
regions,  rsnqinq  ircirn  extreme  semiclassical  to  pure  quantum  ones  in 
order  to  assess  the  modifications  imposed  by  quantum  mechanics  on  the 
classical  chaotic  behaviour 

')  to  Identify  under  what  conditions  these  modifications  take  the  form 
of  a  complete  suppression  of  classical  chaotic  effects. 


T.  Semiclassical  Theoru  of  Electron  Excitation. 


3.1.  Classical  Dynamics  ot  Electron  Excitation 

In  this  section  we  will  develop  the  classical  theory  of  the 
excitation  of  an  hydrogen  atom  in  a  linearly  polarized  monochromatic 
electric  field. 

Here  and  in  the  following  we  will  use  atomic  units,  in  which  the 
Hamiltonian  takes  the  form: 

H=  P^/2  -  1/r  +  ez  cos'a't  (la) 

where  6  and  w  are  the  field  strength  and  frequency  respectively  and 
the  z-coordinate  is  measured  along  the  direction  of  the  external  field. 
The  classical  dynamics  associated  with  (la)  is  conveniently  studied  in 
parabolic  coordinates  since  the  unperturbed  dynamics  is  separable  in 
these  coordinates.  Accordingly,  action-angle  variables  (  n,,  n^,  m,  x,, 
\i,f)  can  be  introduced  /30/,  in  which  the  Hamiltonian  takes  the  torm: 


H=  -  1/2n2  +  €Z  (  n,,  n^,  m,  X,,  X2)  cos  wt  ; 

(lb) 

ri=n,  +  hj  +  |m|. 

Owing  to  axial  symmetry,  m  (which  is  the  z-component  of  the 
angular  momentum)  is  an  integral  of  the  motion;  therefore,  tl) 
describes  an  essentially  2-dimensional  model. 


The  function  z  (n,,  n^,  m,  X,,  X^i  can  be  expanded  in  a  double  Fourier 
senes  in  the  angle  variables  X,,  X2: 


(n,,  n;,,  rn)  e  ‘ 


(ic) 


The  coefficients  z 
given  by/ 1 8/: 


can  be  found  as  shown  in  Appendix  I,  and  are 


I  ’ 

-1,1-1  ■  ■■  *  I  '  \  -.aJ  ^,'. P/^l  I  ' 

-H,  J'^^(Ut(ki  •■  ((  ,  *  l';i)]  ;  tor  r,  ^  k- 

(Id) 

i'  0  for  k^O 

n,-i  =1 

3n  (n,  -  til)  /2  for  k  =  0 

Here  J^.  ere  Bessel  functions  of  the  first  kind  and  J'^,  their  derivates. 
The  dependence  of  z  on  n,,  n^,  rn  is  emPodted  in  the  parameters 
1.1,,  Hj,  which  are  deiineo  Dij 

W I  ,z  "  I ( f'  1  .:•  *  I  I  ^ 

According  to  standard  seinicTassical  approximation  theory  ''28-'', 
give  semiclassical  values  of  dipole  matrix  elements  for 

transitions  n,,i  ->  n',,!  =  n,_i  ♦  k,,i  .  The  element  Zo.o  ,  which  is  just 
the  average  ot  z  over  the  unperturoed  torus  laoeiied  Dy  n,  .  n^  .  m  . 
yields  the  standard  quantum  mechanical  expression  for  the  linear  Stark 
effect. 

If  the  electron  is  initially  in  an  "almost  one-dimensional  '  state,  i.e., 
in  3  state  with  n,>>  no,  n,»m.  then  in  (la-le)  we  can  assume  u,=  !. 
Po  =  0.  In  that  case,  the  dynamics  will  be  described  in  first 
approximation  by  the  one-dimensional  Hamiltonian 

H  -  -  1/2  rr’  +  c  n^’  cos  ■y-r  [3/2  -  2  i  ^  '‘-'■ 

which  is  just  the  Hamiltonian,  in  action-angle  variables,  for  an 
electron  moving  along  the  positive  z  axis/ 1 7,20/; 

H  -  p-  'J  -  l  -'z  *  tz  cos  -yt  .  z  '.I  '2ai 

Wo  start  our  analysis  with  this  simplified  Hamiltonian  (2).  Later  in 
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itii;  lectiufi  we  jhall  discuss,  the  validitii  ol  this  cirie-dimenMunal 
jppro.iimation.  i.e.  we  shjl  discuss  to  what  extent  the  one-din  ensional 
hamiltonian  (2)  is  adequate  in  order  to  describe  the  evolution  ot  quasi 
one-dtmensional  initial  states  under  the  lull  Hamiltonian  (i;. 

Under  appropriate  conditions,  the  classical  system  described  by  the 
Hamiltonian  (2)  undergoes  a  transition  to  chaotic  dynamics.  By  this  we 
mean  that  a  deep  change  occurs  in  the  nature  of  orbits,  which,  above  a 
certain  perturbation  strength,  become  extremely  sensitive  and 
complicated  and  wander  erratically  in  phase  space.  This  irregular 
motion,  if  described  in  the  unperturbed  actions  space,  has  a  diffusive 
character  and  leads  to  fast  ionization.  Quantitative  conditions  for  the 
onset  of  chaotic  dynamics  can  be  obtained  by  rneari'^  of  the  resonance 
overlapping  criterion  / 17,20/.  The  starting  point  of  this  analysis  is 
realizing  that  the  external  field  will  more  effectively  perturb  the 
undisturbed  motion  at  first  order  resonances,  i.e.,  at  values  n  of  the 
unperturbed  actions  such  that  the  external  frequency  co  resonates  with 
some  harmonic  of  the  unperturbed  electron  motion.  These  values  of  n 
are  such  that  si'i  (n)  =  u>  with  s  an  integer  and  i7(n)  the  angular 
frequency  (Kepler  frequency)  of  the  unperturbed  motion: 

Q  (n)  =  dH,)/dn  =  1  /n^ 

First  order  resonances  are  then  given  by  n.  =  (so.i'')'-''^.  However, 
despite  the  the  fact  that  for  these  values  n^  the  perturbation  is  very 
effective,  as  soon  as  it  manages  to  drive  the  motion  away  from  one 
unperturbed  resonant  orbit  its  etfect  becomes  weaker,  ana  nonlinear 
stabilization  may  occur.  In  that  case,  the  motion  keeps  in  a 
neighborhood  (a  "resonance  region")  of  the  original  unperturbed  orbit. 

However,  if  the  perturbation  is  sufficiently  strong,  the  motion  can 
be  driven  so  lar  awau  from  the  original  resonant  value  oi  the  a'  tion, 
that  It  can  fall  under  the  influence  of  another  nearby  resonance.  There, 
the  same  process  may  repeat,  so  that  the  orbit  rn.ay  wander  in  action 
space  in  a  diiiusive  way. 


A  quantitative  estimate  of  the  perturbation  strength  which  is 
necessary  in  order  that  this  may  happen  is  gotten  by  evaluating  the 
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tie.ii'by  regions  overlap  .■■'  I?/. 

riie  -iihilgsis  just  outhnea  can  De  .applied  to  model  i2);  it  is  then 
found  that  for  Vq  -  win'"  >  I  (  where  no  is  the  initial  value  of  the 
action)  and  for  field  strength  exceeding  a  critical  value 

CO  =  trio'*  >  e,.r  -  1/(50  ato  )  (?) 

all  resonance  regions  corresponding  to  n.2;  Oq  do  overlap.  Then  an  orbit 

leaving  with  action  nq  in  a  region  of  phase  space  where  both  fOo  >  I 
and  (3)  are  satisfied  will  diffuse  indefinitelg  and  eventuallg  ionize. 

Notice  that  in  (.3)  we  have  introduced  rescaled  values  c,):  crio"'  for 
Held  and  0.10-  wno"'  tor  irequency.  The  useiulness  ot  this  scaling  is  due 
to  the  fact  that  classical  dynarnics  depends  on  n,-)  only  via  these 
variables  since,  as  can  be  readily  checked,  changing  the  initial  no  by 
some  factor  will  change  the  solution  n(t)  at  any  later  time  by  the 
same  factor,  provided  Co  and  ojo  are  kept  constant,  and  time  is 
measured  in  periods  of  the  field  (see  also  /15/). 

We  emphasize  that  estimate  (3)  is  valid  only  for  (x)o>i.  Indeed  for 
''i.io<i,  i.e.  in  that  phase  space  region  v.'here  to  is  smaller  man  the 
I'epler  frequency,  there  are  no  first-order  resonant  values  of  n,  and 
the  motion  is  therefore  more  stable.  A  transition  to  chaotic  behaviour 
can  still  occur/ 17,20/  due  to  the  finite  width  of  the  resonance  region 
associated  with  mn  =  I  but,  in  order  to  compute  the  in  this  region, 
also  hiqher-order  resonances  Si7  =  pa'  ,  p  >  I  must  be  taken  into 
iccount.  It  is  then  found  that  the  critical  field  increases  with 
dncit'asing  a'li;  however,  for  very  low  static  field  ionization  occurs 
when  -n  ■  .  I 

s'f  course,  higher  order  resonances  play  a  role  in  the  chaotic 
transitiori  also  for  •o,,  1  and,  indeed,  an  approximate  account  of  them 
..  1  already  taten  in  '.  si  via  fiie  choice  01  the  nun.erical  lactor  i/jit. 

'7  A  pi  Opel  reioiid  Ol  der  analysis ■'ZV.-'  leads  to  but  a  small 
iiii  lease  ill  tills,  nunielical  lastni. 
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In  the  chjGtic  regime.  tu'Vr  ’  process  of  diffusive 

excitation  is  convenienttg  described  in  statistical  terms.  Indeed,  an 
equation  of  the  Fokker-Planck  type  can  be  derived  /  1 7/  : 

Of /9  c  =  1/2  9/9n  (D  9f/9n)  (4) 

Where  f  (n.  tl  is  the  distribution  function  and  c  is  the  dimensionless 
time,  measured  as  the  number  of  periods  r='SL't/2Ti  of  the  external 
field.  The  diffusion  coefficient  D  in  quasi-linear  approximation  is 
given  by 

D-  d■'■(.5n)'^>/dT  ^  nViu'y'^^’  nt,)  =  2  nVv^^^  iSi 

t  (4)  and  (5)  were  also  derived  in  /20/).{*) 

'Since  D  increases  with  n  according  to  a  power  la'vv,  it  is  possible  to 
find  an  exact  solution  of  (4).  In  order  to  do  this,  we  must  take  notice 
that  the  stochastic  diffusion  ruled  by  (4)  can  take  place  only  in  that 
part  of  phase  space  where  the  chaotic  transition  has  occurred.  Going 
down  to  lower  and  lower  action  values,  one  will  eventually  meet  an 
invariant  curve  which  has  not  been  destroyed;  we  must  thereiore  look 
for  a  solution  oi  (4)  satisfying  tne  boundary  condition  9f/9n|^^-:  o  of 
zero  flux  across  the  boundary  n-fT of  the  region  of  stability.  In  order  to 
do  that,  the  change  of  variables  y  =  n/Od.  ?  =  TCo-'hi'o'^''^  is 
convenient.  Then,  as  shown  in  app.  11,  for  tV\)<<  I  and  letting  J ^rT /no, 
fhe  solution  assumes  a  suificiently  simple  form: 

ffy,F)-=  (expl-i  l/yy-2/v''^  *  I)’/  F  1  ♦ 

♦  e,'.'pl-(  1/ /y-l)-'/  T  lj/(2u'''‘*/(nF)l  (5a) 


'■•'i  Tfie  numerical  coefficient  2  in  (5'  corresponds  to  the  frequency 
range  1  •-  m,-,  ■  !•.  For  the  asyniptotic  value  of  this  coefficient 

mu’f  be  used,  which  is  near  to  7. 


A  :■  Will  til?  r.Af'tl  III  ieLtlOI'l  ?,  ttll  :■  ICH'IIlUij  iHI'lti  311? ',','lth  th? 

le'jult;.  I.1I  numerical  integration  of  the  equation-;  of  motion  •.■vith 
remarkable  success. 


The  possibility  of  using  this  statistical  description  will  play  an 
important  role  in  our  subsequent  analysis  of  the  ionization  process. 
Indeed,  due  to  the  rapid  growth  with  n  of  the  diffusion  coefficient, 
stochastic  orbits  diffuse  so  fast  towards  high  values  of  n  that  in 
practice  we  can  assume  that  they  actually  ionize-  i.e.  n  becomes 
infinite-  in  a  finite  time.  A  rough  estimate  of  the  ionization  time 
adequate  ror  our  present  purposes  can  be  gotten  from  eq.  (5): 


Cl"  n,j--7D  Cij*-'  )  .  (6) 

In  later  sections  we'll  use  expresslcn  (6)  in  order  to  roughly  estimate 
the  diffusive  ionization  rate  P,  't,"'  • 

In  the  remainder  of  this  section  we  discuss  the  validity  of  the 
one-dimensional  approximation  i2).  Let's  consider  first  the  case  -when 
n,  >>  nv,  n,>>rn  and  therefore  pi<<l.  Then,  since  ~  for 

large  Ikzl,  the  main  contribution  to  the  variation  of  n^  will  be  given 
by  terms  in  (Ic.d)  •/vith  =il.  (Notice  that  give  semiclassical 

matrix  elements  for  transitions  with  An;,  =  ak;.  The  fast  decrease  of 
these  matrix  elements  with  small  U;  when  k;  is  large  was  already 
remarked  in  /3 1/). 

For  ,  the  phase  Xj  begins  to  vary  chaotically,  and  this  leads 

to  a  diffusive  change  in  n-,  also.  The  diffusion  late  for  n-.  in 
quasi-linear  approximation  can  he  derived,  as  -shov.-n  in  /I?/,  by 
retaining  in  fic)  only  terms  with  I;:*!.  One  finds  that: 

0;  ~  ri;  Dill;  ♦  |  m  |  Vti^  i?) 

This  estimate  sho'vvs  that  over  the  ionization  time  ',61  the  change  in 
n-.:  'Ah;!-'  -  ir.  (ri;  ♦  |  tn  1 1  -  n-'  appear-?  to  he  email.  Tiii:  lact 

indicate-?  that  the  onset  of  stocha-sticity  doe-?n  t  le.-jd  to  sigrnticant 


18 


violation;.  01  tlio  ono-dimonsional  apiproximatiori. 


Along  similar  lines,  we  can  show  that  a  suitable  one-dimensional 
approximation  is  valid  also  in  cases  when  n,  ~  |m|»n2.  Indeed  from 

(le)  It  follows  that,  in  such  cases  also,  Then,  upon  neglecting 

Pj  in  (lc,d)  we  obtain  the  one-dimensional  dynamics  for  the  variable 
n,  =  n  -  |m|  ,  described  by  the  Hamiltonian 

H:  -i/2n^  +t  n  cos  (ji)t  [3(n-|m|)/2 

-2U|ri  J'lfIPtk)  k'’  cos  k  X,]  18) 

with  Pt  '5  ( 1- 1 m I /n)'-'^.  v/e  can  now  apply  to  this  one-oimensional 
dynamics  the  resonance  analysis,  just  as  was  done  for  (2).  From  the 
asymptotic  properties  of  J^.'(kP|)  for  k  ->  co  /32/  it  follows  that 

high-k  harmonics  in  (8)  become  exponentially  small  as  soon  as  k  is  so 
high  that  (3/k)^''^  <  m/n  .  This  means  that  the  resonances  of  the  field 
with  such  high  harmonics  (  which  take  place  when  u'o^k  with  (3/k)^''^ 
<  m/n  ),  cannot  significantly  contribute  in  the  chaotization  process. 
Therefore,  when  transition  to  chaotic  motion  is  possible 

only  for  m  2  w'th 

rn^,|.  =  n^i  (3/a'o)^'’  for  a>o>>l.  (9) 

For  (Oq  ~  I  we  may  take  m^^  =5  ng.  At  this  point  we  might  start 

afresh  the  analysis  for  the  Hamiltonian  (6)  in  order  to  determine  the 
critical  field  and  the  diffusion  rate  under  condition  (9).  However  a 
comparison  of  (8)  with  (2)  suggests  that  the  results  of  this  analysis 
shouldn't  deviate  more  than  a  factor  2  from  formulas  (d)  and  (5). 

Again,  the  one-dimensional  approximation  (8)  is  not  significantly 
violated  over  the  ionization  time;  this  can  be  seen  at  once,  because  the 
diffusion  rate  for  O;,  is  still  given  by  (7),  so  that  an  estimate  for  the 
variation  Arii  similar  to  the  previously  established  one  for  (2)  holds 
for  the  present  case. 

Further  details  on  the  classical  dynamics  of  excitation  for  the 


[TiodPl  (2)  Yvill  be  given  in  Sec.  3.2,  where  we  ehall  jlso  dijcij':'.?  the 
results  01  numerical  simulation  of  this  nnodel. 


3.2.  Theorg  of  Quantum  Localization 

The  mam  result  of  the  classical  analysis  carried  out  in  the  previous 
section  was  that  for  sufficiently  strong  field  the  classical  model  (2) 
exhibits  a  transition  to  chaotic  motion.  After  this,  the  classical 
distribution  f(n,  x)  spreads  diffusively  in  action  space,  and  ionization 
fares  place  m  a  finite  time. 

We  will  now  tackle  the  basic  question,  of  what  modifications  would 
be  imposed  on  this  picture  by  quantum  mechanics.  In  particular,  we 
will  study  the  behaviour  of  the  quantum  probability  distribution  over 
the  unperturbed  levels,  which  is  the  quantum  analog  of  f(ri,  t). 

Previous  studies  on  periodically  perturbed  quantum  systems  that 
become  chaotic  in  the  ciaccii^s;  i|„iif.  -  m  particular,  on  the  kicked 
rotator  model  -  brought  into  the  light  the  localization  phenomenon  as  a 
typical  occurrence.  The  quasi-energy  spectrum  is  typically  a  pure  point 
one,  and  quantum  effects  lead  to  a  limitation  of  the  classical  diffusion 
and  to  exponential  localization  of  the  probability  distribution  around 
the  initially  excited  level  np  ;  which  means  that  in  the  average,  and 
apart  from  fluctuations  (that  may  even  be  rather  big  ones)  the 
distribution  looks  like: 

K  exp  (- 2  j  n-rifi  I /I)  (10) 


Here  is  the  time-averaged  population  on  the  unperturbed  level 
rnrrespondmg  to  a  value  n  of  the  quantized  action,  and  1  is  the 
local  I  zation  length. 

In  the  light  of  these  previous  iindings,  H  is  natural  to  assumie  that 
a  similar  picture  applies  also  in  the  present  case,  opecifically,  we 
will  assume  that  even  in  the  serniclassical  region,  and  when  the 
classical  motion  is  chaotic,  a  mechanism  of  quantum  limitation  of  the 


Lhjotic  dilfui-ion  15  worl  inq,  dtid  tiiJt,  under  euitable  conditicine.  this 
mechanism  will  produce  a  situation  analogous  to  the  rotator  case. 
Under  such  conditions,  the  part  oi  the  q.e.  spectrum  relevant  to  our 
analysis  will  he  quasi-discrete;  the  small  line  breadth  of  its  levels 
will  be  negligeable  on  a  time  scale  short  In  cornparison  with  the  very 
long  one  associated  with  multi-photonic  ionization.  While  it  remains 
true  that  the  quantum  atom  described  by  (2)  will  eventually  ionize,  no 
matter  how  small  t,  nevertheless  on  the  time  scale  involved  by  actual 
experiments  the  localization  phenomenon  discussed  here  will  give  it  a 
remarkable  stability  in  contrast  with  the  properties  of  chaotic  motion. 
The  obvious  premise  that  localization  in  hydrogen  atom  is  related  to  a 
unite  time  scale,  should  not  be  forgotten  throughout  this  paper. This 
assumption  will  be  fully  supported  by  the  results  of  our  numerical 
experiments. 

Under  such  assumptions,  we  shall  presently  determine  the 
localization  length  by  the  simple  method  described  in  /3/.  In  this  way, 
we  will  be  able  also  to  determine  the  quantitative  conditions  under 
which  the  localization  picture  actually  applies.  To  this  end,  let  us 
start  with  the  case  of  homogeneous  classical  diffusion,  that  is,  we 
overlook  the  variation  of  D  with  n. 

In  the  semiclassical  regime,  the  evolution  of  a  quantum  state 
initially  coinciding  with  one  unperturbed  eigenstate  n^  will  initially 
follow  to  some  extent  the  classical  development  of  f(n,c).  Therefore, 
over  the  time  scale  in  which  this  semiclassical  approximation  holds, 
the  spread  of  the  wave  packet  over  the  unperturbed  eigenstates  will 
grow  in  time  according  to  An(t)»:(Dt)'''^’. 

However,  the  discrete  character  of  the  quasi-energy  spectrum 
will  prevent  this  diffusive  growth  from  going  on  indefinitely,  as  it 
•would  in  the  classical  case.  The  time  Cj,  after  -which  the  discreteness 

of  the  quasi-energy  spectrum  will  become  manifest  can  be  estimated 
by  '  M,  where  N  is  the  number  of  q.e.  eigenstates  significantly 

excited  by  the  onginal  unperturbed  eigenstate:  indeed  2n/N  is  just 
the  average  spacing  of  q.e.  eigenvalues  significantly  contributing  to 
the  packet  evolution.  Then,  the  number  of  unperturbed  levels  excited  by 
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Mih  vv  ivV  [loCtt-t  jltni'  the  time  l -■  '•  Till:  ineei: 

thjt,  une  unperturbed  level  coritaiiv;.  N  '  An  itj^)  iij.e.  level';  iiid  fhjf. 

one  q.e.  eigenstate  contairrs'  "  An  trj,;i  unperturDed 

levels.  The  latter  number,  however,  is  the  ma.ximum  spread  attainable 
by  the  wave-pacbet,  i.e.,  it  coincides  with  the  localization  length  l. 
There! ore  vve  get  an  equation  lor  tj, : 

sccAnttf))  •-'r  "ro-  I 

■rrhere  '/re  have  introduced  an  undetermined  numerical  lartor  ,  to  be 
found  by  numerical  experiments  /10/.  For  the  rotator  model,  it  was 
round  0  =5  I.  The  same  choice  for  ix  in  the  hydrogen  atorn  case  would 
yield 


I  '-t  0  (iiq)  s  t[,  (11) 

'Where  D(no)  is  given  by  (5):  Dtrio):  Teri"'nfi^/'jL'(i'^''A 

However  this  result  was  obtained  under  the  a'ssurnption  that 
Dr;const.,'','Vhirh  IS  justified  only  in  that  region  '-where  Instead, 

for  1  ■■  rill,  the  dependence  of  D  on  n  may  'sub-stantially  inodify  the 

localization  picture,  and,  if  the  field  strength  exceeds  some  critical 
value.  It  may  even  turn  out  that  localization  is  not  possible  at  all.  (  A 
similar  "delocalization''  phenomenon  was  investigated  and  explained  on 
a  simple  example  in  Refs.  /3,t0/). 

In  order  to  clarify  how  delocalization  occurs,  we  need  to  modify  tne 
above  method  for  determining  I  ,  in  such  a  'way  that  the  dependence  of 
[i  on  n  I'S  explicitly  taken  into  account.Therelore.  in  place  ot  An('i  c 
'['  I'rii,'  r'l’  we  mu'St  substitute  the  dependence  of  An  oii  t  that  is 
entorreo  by  the  Fol  ker-Planck  eguation  (T).  In  this  'vvay  'vve  find,  as  a 
result  of  the  calculations  developed  in  Appendix  111,  that  An(i:)  is 

ijl'.'eii  fjij: 


COiiri;  ji,  1 


(1 :) 
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By  tlif  ijint'  jrquttit'nf  o;.  ■itiovt'.  yvi?  ljii  now  litid  rj-|  ond  I  iruni  tho 
iod jli33tion  Londition  Atv.'C|^)  =  .  Howevei',  n  C|j  id  large  enough,  the 

curve  An(r)  will  never  intersect  the  straight  line  ar  heiore  explnding 
at  t  =  (jLio"''V3eo^  .  When  this  happens,  no  localization  is  possible  and 
this  implies  unbounded  diffusion  for  the  electron.  More  specifically,  in 
App.  Ill  we  show  that  the  solution  of  An(tp)  gives  the 

localization  length  I  in  the  form: 


1  -4  «■'(  u  /3  (13) 

where  u  is  the  least  of  the  two  solutions  of  the  equation  (*) 


3ot:‘(<jj(|'^'°  €(,)“'  =  g(u)  =  u(  1 -u)‘/(2-u) 


(H) 


such  that  0<u<i.  Numerical  data  indicate  that  here  too,  like  in  the 
rotator  case,  a-'--  1  is  to  be  chosen  (  see  e.  g.  fig.  10  and  related 
comments  in  sec.  3).  Therefore,  since  the  function  g(u)  in  the  interval 
(0,1),  has  a  maximum  ^  1/12  at  u  =  (3-V'5)/2.  it  follows  that  (or 


'0 


a)o^^V/(6no) 


(15) 


eg.  (14)  has  no  solution. 

Thus  defines  the  threshold  for  quantum  delocalization.  Of 
course,  in  order  that  delocalization  may  occur,  it  is  also  necessary 
that  t'n  exceeds  the  threshold  for  classical  chaos  (3),  just  because 
the  semiclassical  estimate  (12)  holds  under  the  assumption  that 
chaotic  diffusion  takes  place  in  the  classical  system.  According  to  the 
argument  just  outlined,  across  the  threshold  i  qualitative  change 


‘  j  The  s  light  difference  in  numerical  coefficients  between  '.loi  and 
the  analognijs.  lorrnula  oi  rei.  //I/  is  due  to  the  lact  that  in  ,'j|;  j 
value  of  0  soriiewnat  le-.ser  than  i  vvas  introduced  m 


M. 


uccui'i,.  onJ  the  locdli^aticifi  picture  is  no  lonqer  lUc-titied.  We  shoulO 
then  expect  that  above  this  threshold  here  is  no  quantum  limitation  to 
the  classical  diifusKm  and,  indeed,  this  will  clearly  appear'  from 
numerical  results. 


The  above  one-dimensional  analysis  can  now  be  modified,  so  as  to 
apply  also  in  the  2-dim  case  for  quasi- 1-dim.  states.  Indeed,  even 
though  vve  know  from  Sec.  2.1  that  the  1-dim.  approximation  is  not 
violated  for  such  states,  still  we  cannot  i?  ;t7/'7t7/7  exclude  that  the 
presence  of  an  additional  degree  of  freedom  can  destroy  the 
localization  of  these  quasi,  but  not  strictly,  1-dim.  states.  However, 
we  can  answer  this  question  by  the  same  method  used  in  the  i-oim. 
case.  Indeed  since  classical  diffusion  will  now  take  place  both  for  n, 
and  n^  ,  the  number  of  unperturbed  levels  excited  at  time  c  will  be 

M  ~  A  n,  (r)  A  nj(t) 

For  Arii(t)  we  will  now  take  eq.  (12);  moreover,  since  An^fT)  can  be 
assumed  to  be  small  '.see  sec.  2.1)  ,  we  will  take  AnjU)  "  (Dvi:)''^' 

with  D'j  as  in  formula  (?).  Imposing  now  the  delocalization  condition 
N  1  r;  /3/6/10/1  i/ifi/2i/  one  easily  get  the  estimate  for  the 
two-dimensional  delocalization  threshold: 

e’o  '■  ^  [n2  (n^  *  |m|)l'^‘  (16  ) 

where  again  an  undetermined  numerical  factor  3"'  was  introduced.* 

The  estimate  (16  )  clearly  Indicates  that  2-dimensionality  sharply 
decreases  the  delocalization  threshold.  Nonetheless,  for  states  with  m 
-n^  -■  1  the  2-dirn.  threshold  is  almost  the  same  as  the  1-dim.  one. 


(*)Actually  a  more  refined  analysis  shows  that  in  the  2-aim.  case  here 
considered  one  has  localization  on  an  exponentially  large  scale  i6ee 
ref.  qT  and  eq.  (T.^i  of  ref.  61.  so  that  delocalization  takes  place  only 


j-i 


V/e  jce  therefore  justified  m  assuining  that  for  such  quasi- 1 -dim. 
states,  the  localization-delocalization  picture  remains  valid. 

The  decrease  of  the  delocalization  threshold  in  the  2-dim.  model 
yields  one  reason  (another  one  will  be  given  in  Sec. 3. 3  )  why  some 
agreement  was  found  between  experiments  on  excitation  of  H  atoms 
from  states  with  n~66  with  field  frequency  w/2‘n''  10  GHz  /12/  and 

the  results  of  numerical  simulation  of  this  process  on  the  c/jss/cj/ 
model  /1 5/26/.  Indeed,  in  numerical  experiments  the  initial 
distribution  of  states  with  n-b6  was  nearly  microcanonical,  so  that 
the  above  discussed  2-dim  effects  played  an  essential  role  in  lowering 
the  delocalization  threshold. 

(n  closing  this  section,  a  couple  of  remarks  conrerning  the  validity 
of  the  quantum  system  described  by  the  Hamiltonian  (2)  as  a 
physically  realistic  model  are  in  order.  In  the  first  place,  in  the 
quantum  theory  discussed  above  we  considered  the  electric  field  as 
classical.  This  approximation  holds  if  the  full  number  N  of  field  quanta 
inside  the  microwave  cavity  of  volume  V  : 

N  :  /(4Ti1iio)  ~  3.10‘’€o^''^ho'''J‘'o) 

IS  sufficiently  large.  For  istance,  for  e,,  :.05,  U'o  =1,  n^  :I00, 
V:icm\  which  are  typical  for  the  range  explored  in  our  investigations, 
we  get  N  ■'ID". 


Also  the  question  may  be  raised  whether  the  diffusive  excitation 
process,  that  Is  made  possible  by  the  delocalization  phenomenon, 
should  not  be  significantly  reduced  by  the  spontaneous  radiation 
process.  However,  the  rate  f,,  of  the  latter  process  is  much  less  than 

the  diffusion  rate  r^,.  indeed  even  tor  orbital  quantum  number  l  i  the 
rate  P,  "(c’n,-,''  1^)"'  i  c  /3z/.  Estimating  P^,  by  the  inverse  oi 
the  classical  ionization  time  i.e.,  by  r'|''ci'/2ii.  with  r|  a.-,  in  (6),  we 
obtain; 


r/r^,  -  iM-^ 

Where  the  numencal  estimate  is  given  tor  the  tgpical  values  Wfi-l, 
eo-.OS;  notice  also  that  the  ratio  does  not  depend  on  Oq.  Actually,  this 
ratio  is  even  smaller  because  the  extended  state  contains  1  up  to  -  /ri,j 

»l. 

3.3.  Ionization  in  the  Presence  of  Localization 

According  to  the  theory  developed  in  the  previous  section,  as  long 
as  the  one-dimensional  approximation  is  valid,  the  dependence  of 
ionization  on  the  field  strength  should  have  a  more  or  less  marked 
threshold  character,  denned  by  the  quantum  delocalization  border  (15). 
However,  a  microcanonical  distribution  of  initial  states  looks  fairly 
typical  in  many  physical  situations,  so  that  it  is  interesting  to 
investigate  what  should  be  in  that  case  the  dependence  of  ionization 
probability  on  the  field  intensity.  Indeed,  since  the  tvvo-dim. 
localization  border  (16)  depends  on  both  quantum  numbers  n,.  rij,  in 
that  case  we  should  expect  that  for  any  (  not  too  "high”)  field,  a 
traction  of  the  states,  depending  on  the  field  strength,  will  be 

delocalized,  while  others  will  be  localized  and  will  therefore  give  ho 
contribution  to  the  ionization  rate. 

We  will  derive  this  dependence,  under  the  assumption  that  the 

interaction  time  of  the  atom  with  the  field  is  large  enough  for 

the  classical  system  to  undergo  complete  ionization  i.e.,  that  the 

classical  ionization  probability  Pj’^'^l.  Besides  that,  however, 

must  not  be  so  large  that  direct  quantum  ionization  from  the 
stationary  distribution  (lO)  into  the  continuous  spectrum  becomes 
effecti'-'etsee  the  comments  in  sec.  2.2). 

Let's  first  assume  that  we  have  initially  a  homugeneous  distribution 
of  states  with  a  fixed  value  of  Of,  and  of  the  magnetic  quantum  number 
III.  Then,  after  the  time  all  atoms  initially  in  states  with 

rii>  !f^'*\iv.,m)  will  be  ionized. These  are  precisely  tiie  atoms  initially 
in  states  with  n,  >  w’th  n/'"  given  by  the  equation 


26 

Then,  recallinq  that  rii+n^  =  ni;i-|m|,  we  see  that  the  fraction 
of  atoms  m  the  ensemble  which  will  not  be  ionized  at  time  is 
equal  to  I'n,'^'"  /(oq  -|m|)  (the  factor  2  is  due  to  symmetry  for 
exchanges  n,  02).  Computing  n2''’  from  eq.  (16),  we  get: 

=  [(m"  +no‘A')' -  |  m  |  ]  /2  (17) 

The  ionization  probability  is  therefore  given  by: 

Pp  1  -  2n2'^7(no- 1  m  | )  = 

=  I  -f  [m‘+no‘A‘l'''^'|m|}/(no-|m|) 

where 

A=  •  ns) 

Now  let's  assume  that  the  initial  distribution  is  microcanonical, 
i.e.,  that  all  quantum  states  with  a  fixed  0,5  are  equally  represented  in 
it.  The  full  number  of  such  states  is  no^  .  For  any  given  value  of  m. 
the  number  of  non-ionized  states  at  time  for  the  given  Cq  is  Just 

(m,e'o).  We  must  then  sum  over  the  different  values  of  m;  in  doing 
this,  however,  we  must  remember  that  there  is  a  classical  value  rn*^’’ 
above  which  there  is  no  ionization  (9).  Replacing  the  sum  by  an 
integral,  we  find  that  the  fraction  of  non-tonized  atoms  at  time 

for  the  given  eg  is  given  by: 

1  -  Pj  %  (4/no^)/  02'^''  (m,  Cfi)  dm 

The  factor  4  in  the  above  formula  is  due  to  symmetry  'with  respect 
to  exchanges  m  -m,  n,  n,  .  The  latter  symmetry  must  be  taken 
into  account  also  in  inserting  the  appropriate  expression  for  02’''^  in  the 
integrand.  Iride'‘d,  cannot  e:;cecd  (n,-)- |  m  |  )/2:  otherwise,  since  the 
argument  is  symmetric  in  n|,  rij,  a  supercritical  value  of  02  would 
enforce  a  subcntical  value  of  npno-|m|-n2.  Therefore,  n®''  in  the 
above  integral  is  actually  the  intimum  between  (17)  and  (no-|m|)/2, 
i.e.,  it  IS  given  by  (17)  for  |m|<mi  =ri5v'(  1 -A-),  while,  for  m>m,  it  is 
equal  to  (n- 1  m  |  )/2. 


TtiPn.  dssuminq  -3  Oy  (which,  as  we  have  already  remarked,  r? 
legitimate  for  uyj  1),  and  evaluating  the  integral,  we  finally  get  the 
dependence  of  Pj  on  the  field  Cy  in  the  following  form: 

P,  :  /(1-A^)  -  ln{  [I+/(1-A2)1/A^  }  (19) 

where  A  is  given  tiy  eg.  ( 18). 

Unfortunately,  it  would  not  be  correct  to  use  available  experimental 
data  as  a  check  of  (19),  for  the  following  reasons. 

In  the  first  place,  whereas  experimental  data  concern  the  frequency 
region  (On<1,  the  above  described  theory  of  localization  was  derived 
in  the  frequency  region  (jL)y>1.  where  1st  order  resonances  exist  (the 
peculiarities  of  the  excitation  process  for  Wy  <1  will  be  discussed  in 
sec.  5.2). 

Second,  according  to  numerical  data  /1 5/,  in  experiments  the 
condition  P,'='=l  was  not  fulfilled  after  time  ;  indeed,  by 

increasing  a  further  increase  of  P,"^'  was  gotten.  This  fact 

makes  impossible  the  comparison  of  available  experimental  data  with 
(19). 

3.4.  Comparison  of  Diffusive  and  One-Photon  Ionization 

In  the  delocalization  region  the  quantum  mechanism  of 

suppression  of  classical  diffusion  is  not  at  work  and  therefore  one 
expects  that  the  quantum  electron  will  diffuse  and  ionize  like  the 
classical  one.  This  fact  has  been  numerically  checked  and  will  be 
discussed  in  section  3.  The  resulting  diffusive  excitation  can  hardly  be 
described  vnthin  the  framework  of  conventional  multiphoton  theory; 
moreover,  it  usually  takes  place  in  a  very  different  range  of 
frequencies  than  considered  there.  In  order  to  appreciate  the 
etierriveness  oi  rtns  new  ionization  process,  t  is  interesting  to 
compare  it  with  the  familiar  one-photon  process. 


T Cl  obtain  a  quantitative  estimate  lor  une-[.inotun  ionisation,  we 
shall  first  observe  that  any  normalized  energy  eigenfunction  for  the 
unperturbed  one-dimensional  hydrogen  atom  (i.e.  tor  the  Hamiltonian 
(2)  with  c=0)  can  be  wntten  as  u(z)  =  z  R(z),  where  R(z)  is  a  radial 
eigenfunction  for  the  3-dim.  atom  with  orbital  quantum  number  1=0. 
Therefore,  the  matrix  element  for  the  photoelectric  transition  from 
the  n-th  unperturbed  level  of  the  I -dim.  model  (2)  to  the  continuum 
state  having  energy  P^/2  =  -l/2n^  »  is  given  by 

--  |dzz^R;/z)R^(z)  (20) 

where  Rp^  are  radial  eigenfunctions  for  the  3-dimensional  atom 
with  orbital  quantum  number  -0  (  We  assume  Rp,^  to  be  normalized  on 
the  energy  scale). 

For  highly  excited  states  n>>l  the  integral  (20)  can  he  evaluated  by 
semiclassical  methods.  In  Ref.  /33/  the  following  semiclassical  value 
01  dipole  matrix  elements  for  transitions  from  states  (n.i)  to  (p,  l+i) 
was  found  for  l«n,  p«1  (notice  the  difference  in  normalization 
between  (20)  and  ref..''33/): 

R^'’'‘'n,i  1  K,..j  (CO  lV3)l  /(3 ncon’''^) 

where  (^)  are  Mac  Donald  functions.  Considering  that  the  2nd  term 
in  square  brackets  is  negligible  for  small  1,  and  that  for  £,  ->  0, 

s:  0.459(3n')''’^'(3^/2)'^'^  we  find  the  follov,'irig  semiclassical 
value  for  (20): 


RP,0  Sj  p;P.I  S;  0.451:} 
ri,o  riu 


(-1)  n  (j.r 


Then  the  transition  probability  for  unit  time  is: 


r,j,  n/2)  ^  'J"  t.5  c-/(>o"^'’  'I 


and  i.he  ionization  probability  in  one  period  of  the  external  Held  is 


1 


(tor  .'j.Ti  >  rifl/j).  This  value  is  ^-/j  times  'arqer  than  in  ret.  /''T/  Oue  to 
averaqinq  over  solid  jtiqle. 

In  order  to  compare  the  l -photon  ionization  and  toe  duiusive 
ionization  we  shall  choose  tne  optimal  regime  oi  each  process.  Then 
for  1 -photon  ionization  we  take  <x)~(2n, so  thatT|,=r  34  t|j-’ 

Instead,  'or  diffusive  ionization  we  tai. e  o.i  n„  ’  "  i).  and  vve 

estimate  ttie  ionization  probaPility  per  period  as  T^,  s  t]''with  tj  as 
in  ioi: 


In  this  wau  we  see  that  diffusive  ionization,  wiiicfi  fates  place  for  a 
much  lower  frequency  than  I -photon  ionization,  is  a  much  more 
eifective  process  than  the  latter: 


In  real  physical  time  this  ratio  changes  as  each  i  is  multiplied  Dy 
Its  oven  optimal  treguency  which  gives: 

rpcr.j,  -n./'Vs  (24) 

This  ratio  is  still  large  for  nn>'^  i- 

A  detailed  analysis  of  the  dependence  of  the  ionization  proPahility 
on  frequency  will  he  discussed  in  sec.  3.4 


'.-I  1 1  ih  I j  f  1 1  ill  tiij  7  milit-l  lliQ  jiij  t  flijijt'h  p  jldliiii't  el . 


we  snail  now  oisruss  some  peculiarities  ot  fhe  exciration  ana 
ionisation  for  quasi-classical  states  in  the  classically  stable  region 
a)|j  >  1.  In  this  case  the  classical  motion  from  one  resonance  to 
the  next  one  is  lorbidden  by  the  presence  oi  smooth  invariant  curves 
between  them. 

Due  to  this  fact,  excitation  and  ionization  can  take  place  only  thanks 
to  tunneling  through  the  classically  forbidden  region.  A  distinction  can 
then  be  made  between  the  two  opposite  cases,  when  the  number  of  levels 
coupled  by  the  field  in  one  period.  An  D'-'*  "tii  ni^/'i'i-i"”''  is  small 

or  large.  In  the  former  case,  perturbation  theory  holds.  Therefore,  the 
probability  of  transition  from  the  initially  excited  level  no''  I  to  nearby 
levels  IS  small,  and  the  ionization  probability  depends  aloobraically  on 
Co- 


w,  o  (no  Co/'Oo''^')"'^'  (25) 

where  k  '  no  is  the  number  of  photons  required  for  the  transition  into 

the  chaotic  component.  In  the  opposite  case  when  An>>l  but  still 
tunneling  becomes  domimant.  and  we  can  reasonably  expect  Wj  to 

depend  on  Cy  according  to 

Vv'j  i>  '|  *  expl'  cno  (Cj.^  ■  (26) 

.vtiere  i.  is  3  numerical  constant  of  the  erder  of  unitu,  r|‘'represents  the 
lonizaf iiin  rate  irom  me  cnacifir  rcimponent,  and  fhe  e:-:poneritiai  factor  is 
"1  i-suKieii  111  funupiitig  mfn  me  ci.assically  lorbidden  region  ,  on 

account  ui  It:  analogy  with  the  formula  describing  tunneling  in  a  static 
lit'M.  ■’’tiM  iirnlitiiiM  frir  ippl ii‘ atu 1 1  ty  of  (26)  is 

,  -  l/iSUvo,/"’  )  (27) 

111  ui  del  mjf  nuhJ,.,  li'ViM..  ii|,  '  Di'-i'y-''-  ai  e  regull  ed;  e.g.,  lui 


r 


■  1 


'I'n  I  If'.vi:  with  hi,  '  5''  IIIN:t  Ll'>  I.  Ull  :lO''l  f'lj. 

it  I'i  interesting  to  compare  these  results  with  t.eldgsh's  theory  lor 
tiinneiing/7-1/  in  which  an  adiaoatic  parameter  T  =  ni/tcnu)  =  'c.ij/e,-, 
IS  ititroduced.  discriminating  the  perturbative  regime  ':Y>>1)  from  the 
tunneling  regime  (Y<<1).  In  the  present  case,  (27)  shoves  that  in  order 
that  tunneling  loniration  according  to  (26)  can  take  place  it  ic-  necessary 
that  Y  -  a'o/'cn  >>1  while  for  Y<<l  one  falls  into  the  region  where 
diffusive  ionisation  occurs.  For  iristance,  tor  uh-i  ■■  I,  Cn''  0,01,  rii., 
'''•100  eg.  (2?)  is  satisfied  but  T  s:  lOO.  Therefore,  we  see  that  here 
also,  like  in  /34/,  the  rnuitiphoton  regime  occurs  for  weak  field  in  the 
perturbative  region  (>;„  ■^^<p).  whereas  tunneling  takes  place  in  tne 

uppusite  case  of  strong  field  (c,-,  •'’■kp).  In  i;unc lusion,  due  to  the 

phenomenon  of  dilfusive  excitation,  the  l'eldi':ii  parameter  loses  it; 
usual  meaning  and  a  new  parameter  Y^  must  be  introduced  in  order  to 

discriminate  between  tne  perturbative  regime 1  land  tne  tunneling 
regime  (  rp|  -r|  ).  Accordinq  to  the  previou-s  discussion  the  new  parameter 
Y^  will  have  the  expre-ssion 


(hfiCn)  - 


•  Pro 


/(rn, 


(26) 


4.  Numerical  Results 


4.1.  Methods  of  Numerical  Simulation 

In  this  section  we  shall  describe  the  numerical  methods  and  the 
checking  procedures  -we  used  in  our  computer  simulation  of  the 
classical  and  quantum  dgnarnics  of  the  l-dirn.  model. 

Reducing  to  one  the  dimension  of  the  problem  sharply  decreases  the 
computation  time  in  the  quantum  case,  and  this  allows  tor  a  more 
precise  investigation  of  the  e.xcitation  dynamics. 

The  main  computations  were  carried  out  on  the  CRAY-XMP  Computer. 

The  numerical  solution  of  the  classical  equations  was  carried  out  in 
action-angle  variables  (n.  X)  .As  in  /1 8/,  in  order  to  circumvent  the 
singularity  at  2=0  a  change  was  made  to  new  variables  (n.  t  )and  to  a 
new  time  T|  ,  which  allowed  to  write  the  equations  in  the  following 
torin: 


dn/dT|  --  -  f.  r\^  cos  out  sin  ^ 

d.^/dt  :  n"'  ♦  I’en  cos  out  (1-  cos  0  (29) 

dt/dT|  =  1  -  cos  C,  ;  X=  ^  -  sin  c, 

M  similar  method  for  avoiding  the  singularity  at  the  origin  was  used 
in  /  29/.  Eqs.  (29)  were  then  numerically  integrated  by  the  Runge-hutta 
mefhod.  The  initial  distribution  of  classical  trajectories  was  taken  on 
a  line  in  phase  space  with  n=n(,  and  uniformly  distributed  phases  X; 
this  choice  corresponds  to  the  initial  condition  used  in  the  quantum 
case  (only  one  level  excited  wifli  n=nn  ).  The  full  number  of  classical 
traiec tones  vvas  t.aien  jv'j  oi  I'i""!'. 

^^n  absorption  mechanism  was  introduced  for  trajectories  being 
excited  atiovc  n  "r  dnii.  A  ifi.jnqe  in  tne  border  oi  absorption  only 
vveally  aiiei  ted  the  evcifafion  probability. 

The  investigation  of  guanturn  dynamics  described  b'j  the  Hamiltonian 
(2)  was  carried  out  by  t,vo  distinct  methods.  In  the  first  one. 
followinq  /  1 8/,  a  base  of  discrete  unperturbed  eigenstates  was  used, 
and  the  equations  were  solved  for  the  amplitudes  c^  of  the  expansion 


•j I  ••tott'  vH!*  tot  OVOt'  thOsK  OlOt’fri'tdt Oi'J 

ir-  =  -(t/jn^)  c,  +  t(t)  7  (30) 

f\  n  —  fifi  fi 

Tlie  Vdlue  of  was  approxlmatelu  20-40  lovols  lesser  than  the 
initially  excited  state  n^  .  A  further  decrease  of  did  not 

appreciably  influence  the  dynamics,  owing  to  the  exponential  decrease 
01  the  distribution  -  |Cf,|^  in  the  region  n<no  where  the  classical 
inotion  is  stable.  A  typical  value  for  the  full  number  of  levels  for 
which  eqs.  (30)  were  solved  was  ND=  n^^._.  -  n^,^  =192. 

In  order  to  numerically  integrate  (30)  the  time  dependence  of  the 
field  was  approximated  by  c(t)  =  At  e  cos  u>t  7^^.  6(t  -  l  At)  with 
At  =  2ti7ojL,  where  L  is  the  number  of  integration  steps  per  period. 
This  scheme  of  integration  is  physically  equivalent  to  introducing 
supplementary  fields  with  frequencies  =  kLo.!,  k=l,2,  ...  Since  in 

our  computations  ai"  l/rio^  and  the  number  L  of  steps  was  chosen 

between  lOO  and  500,  then  even  the  frequency  cO]  ^  lOOco  v/as  much 
larger  than  all  frequencies  for  transitions  between  intermediate 
levels.  Therefore,  the  influence  of  the  fictitious  frequencies  o.>^  can  be 

considered  to  be  small. 

The  integration  of  the  numerical  scheme  thus  obtained  can  be 
carried  out  exactly;  indeed,  It  reduces  to  successive  applications  of  a 
matrix  to  a  vector  c(t)'. 

c(t  +  At)  =  T  exp[  -  icn  (cos  a>t^)At  z]  c(t) 

=  T  Q  Z  Q''  c(t) 

where  T  and  Z  are  unitary  diagonal  matrice':.,  with  T^^  =  expiiAt/fxn") 
and  Z|,|  ^  =  expl-Uo  (cos  a'tj,)At  zj,  z^^  are  the  eigenvalues  of  the 
iTiatrix  z^i  and  Q  is  a  unitary  matrix  that  carries  the  same  matrix 
into  diagonal  form.  Vv'ith  this  procedure,  tfie  normalization  Vv'=j  Unh 
=  1  IS  conserved  to  a  very  high  accuracy  r  icC'  i.  In  / 16/21/  the 
iperator  e:<p(-i AtCuZcosa'^j  was  computed,  by  means  of  its  expansion 
in  powers  of  At  (up  to  the  5th  order),  which  led  to  an  effective 
damping  on  higher  levels  and  to  a  poorer  conservation  of 


3^^ 

iiOlTfijl  IZit lOti.  The  rieYi'  inethoil  H'r-ed  here  .jppejr'.-.  -.-.iciriU  tc  jritly  Inure 
el I'lL'ieiit.  in  that  it  permits  to  decrease  the  niimtier  of  steps  per 
period. 

The  main  inconvemente  with  the  just  described  integration  scheme 
is  that  the  continuous  spectrum  is  completely  neglected.  Even  though  a 
number  of  arguments  can  be  put  forth  /1 6/21/,  suggesting  that  the 
continuous  spectrum  would  not  essentially  modify  the  dynamics  of 
excitation  over  discrete  levels,  nevertheless  it  is  important  to  build  a 
numerical  model  free  of  this  shortcoming. 

As  tar  as  we  know,  no  numerical  experiments  vvote  u[i  to  novY 
performed,  giving  a  precise  account  for  continuous  spectrLirn/S  1 /.  A 
partial  consideration  of  transitions  into  the  continuum,  has  been  given 
in  ref.  /36/.  However  no  account  was  there  taken  for 
continiinm-continuum  transitions,  which,  generally  speaking,  do  not 
appear  negligible  as  compared  with  transitions  to  and  from  the 
continuum.  Moreover,  the  number  of  equations  to  be  solved  sharply 
increases  with  the  level  number  nj  and  this  does  not  allow  for 
investigation  oi  excited  states  with  no  60. 

A  more  efficient  account  for  continuum  can  be  given  by  means  of 
the  so  called  Sturm  base.  This  base  is  introduced  by  considering  the 
following  eigenvalue  equation: 

-  1/2  d^u/dz^  -  (^/z)  u  =Eu  z>0.  E<0,  ^>0  .  (31) 

For  J.=  l,  (31)  IS  just  the  Schrodinger  equation  for  the  stationary 
states  of  the  unperturbed  I -dim.  hydrogen  atom  .  By  changing  variables 
according  to  ^=2z,  u(z)-  (^/2)'''^v  (0,  eq.  (31)  becomes 

SjV  =  d/dC  (F  dv/dO  ♦  |(E/2)  c  -  1/(4.-;)  Iv  -  -fw  (32) 

The  Sturm  base  is  generated,  by  considering  (32)  as  defining 
eigenvalues  -fl  for  the  operator  Sp,  where  E<0  is  an  arbitrary  fixed 

parameter. 

Instead,  considering  in  (32)  f.  as  parameter  and  E  as  the 
eigenvalues,  one  would  resort  to  the  usual  base,  including  continuum 
eigenfunctions. 


'S 


it.  i.  Inuvvti  .'JS/  thJt  S(^  IS  j  splf-jdinint  opprjtnr  vYith  j  puiely 

discrptH  ipHCtrum  f-,.=  (s*  l)(-2E)''",  with  s>ij  jn  inteqpr. 

EigHfiturict ions  tor  5^  ore  given  Dy 


F(-s.  2.  £(-2E)''-i  li(-2E}'-l''- 
exp|-(2/2)t-2E)'''‘l 


iiid  ore  orthonortnol: 


Here,  jiid  hetovv,  f  '.vith  three  vdrijCdee  ■.■'.ill  inditiie  tiie  otniiluent 
hypergeornetric  iiinction.  In  the  following,  we  ‘rii.ill  chuose  [--1/20,1', 
iVi  heinq  the  initially  excited  level.  Then,  i^^^_,t2ri  ic,  apart  from  a 
normalization  constant,  the  On-th  unperturhed  eigenfunction. 

We  al’so  need  matrix  elements  for  2  and  c*‘.  For  2'-  they  are  given  in 
/28/,  vvhere  they  are  used  in  order  to  calculate  the  2nd  order  Stark 
effect,  rtatrix  element's  for  .?  can  he  ohtair.ed  hy  direct  cornputation. 

Non  zero  elements  for  2  and  are  then  given  hy: 

^5,5  ■  2nn(s*  1 ) 

t.i' Drin^  (S*  1 )/ 

-2nn/(2sH  il'Sc 

'2-':y  ^  '2-i^  T  n,,-s./'S'^-  I; 

Let  nu'(Y  't'tf i=  ■.•'z  't’(t)  he  the  solution  of  the  iSchroedinger  equation 

.  V 1 1  h  the  H.jiiii  1 1  iiiiian  '!2i,  and  .''-2  '2.  Then,  l-'H  [  '  |'i'|'d.2-t.  Since 

*  I 

the  I.  ;  mat  e  up  a  complete  oi  thonorinal  set,  t  can  he  expanded  m  the 


form 


with  Eo=  -  l/2rii3‘  and  'HoM  / 2/r\(,  j corresponding  to  the  irntiany 
excited  level  rif,.  By  using  the  orthonorrnahty  of  the  f^  and  the 
expression  (33)  for  matrix  elements,  from  the  Schroedinger  equation 
we  obtain  equations  for  the  amplitudes  A.(t): 

2(sH)A5  -  [s(s*l))''"4,_,  -  l(s*IMs+2)l''''^  A^^P 

:-i(l2(s+ l-no)/n(|2l  a^  3e(t)no(s+ 1 'r  A^^  (34) 

-  e(t)riQl(2s+ l)(s(st  1))'''^  A^_l  *■  (2s+3)((s+ l)(s-*2))''^-^ 

+  (ri,/2)e(t)[s/($^-U  A,^.2  ♦  (s-2)  ((s+2)2- l)i’'2  a,.,!) 

This  infinite  system  of  equations  is  exact  and.  even  though  only  a 
discrete  base  was  used,  it  completely  takes  into  account  the 
continuum.  Indeed,  each  Sturm  function  is  a  superposition  of  several 
eigenfunctions  from  the  unperturbed  base,  including  eigenfunctions 
belonging  to  the  continuous  spectrum. 

Once  eqs.  (34)  have  been  solved  for  A^,  the  original  amplitudes  c.,|(t) 
of  the  expansion  of  fit)  over  the  unperturbed  eigenstates 

u^(r)=  2zir^''^  e'  F(-n+l,  2,  2z/n) 
can  be  recovered  by 


'.'.here  fhe  it  jncformation  matrix  B,.  irom  the  Sturm  to  tho 

rit. 

unperturbed  base  is  given  by 


:  I/.'  I  d/(t;'2'''‘’t T-'i  iy2/2)“ 

f  f 

;4(:-.*  1  wi|)i  (zio*  Dili'  '*’  (M|,-n)'*'’  *  (nijn/-  I  i'‘  l(iy,*ni' I 
F(-s,  -(irl),  2,  -  drifin/fn-rin)') 


Here  F  is  Gauss'  hypergeiametric  function.  A  similar  computation,  for 
continuous  spectrum  unperturbed  eigenfunctions  can  be  made,  by 
simply  substituting  i/p  in  place  of  n,  p  being  the  electron  momentum 
(an  analogous  method  was  used,  e.g.,  in  /33/X  A  method  for  Lornpucinq 
F  with  large  s,  n  is  given  in  3ppendi.x  IV. 

The  numerical  integration  of  eqs.  (35)  was  performed  as  follows. 

One  level  Sn-rin-l  was  initially  excited,  so  that  A^(0)-t/2/nf|) 

Then  eqs.  (34)  were  solved  for  s  x  As  a  rule,  1 0-30, 

and  the  full  number  NS  of  Sturm  levels  ranged  from  256  to  576.  The 
dependence  of  the  field  on  time  was  taken  in  the  same  way  as  in  the 
previously  described  method,  with  approximately  the  same  number  of 
steps  per  period:  100<L<500.  Just  as  in  the  1st  method,  the 
introduction  of  delta  functions  into  the  numerical  scheme  made  it 
possible  to  exactly  integrate  the  truncated  set  of  equations  (34)  by 
repeated  applications  of  matrices.  For  the  same  reason,  the  loss  of 
normalization  was  very  small  (  ~  lO"’).  Unlike  the  1st  method,  here  the 

presence  of  high  frequencies  (0^.=  kLo,'  led  to  direct  transitions  into  the 

continuum;  however,  for  the  chosen  values  of  L  the  probability  of  such 
transitions  was  negligibly  small.  For  instance,  for  no=60,  u>d=1, 

L=100  we  get  <jO,no^  =  50  and  'T^  553.10'^  to  be  compared  with 

-  .02.  Therefore  the  small  5-function  kicks  introduced  by  the 
numerical  simulation  of  the  monochromatic  perturbation  do  not  have 
any  effect  on  the  physics  of  the  problem;  moreover,  their  influence  can 
be  kept  under  control  by  varying  the  integration  step. 

In  our  opinion,  monochromaticity  of  the  perturbation  is  important 
lor  this  problem,  and  substituting  a  6-like  perturbation  (t-2Tik/o.)) 
m  place  of  t  cosmt  /37/36/  can  lead  to  a  significant  modification  of 
the  physical  picture  of  rnultiphotonic  excitation.  The  role  of 
multiphoton  transition  in  the  2'-dirn  model  with  a  6-like  perturbation 
was  studied  in  ref.  /57/. 

The  values  of  A^(t)  obtained  by  integrating  (34)  were  used  to  find 
the  amplitudes  c^|(r)  over  the  unperturbed  discrete  base  by  means  of 
the  transformation  matrix  In  this  way  c^(')  were  found  for 
approximately  200  levels.  Since  the  total  probability  was  conserved 
with  high  accuracy,  it  was  then  possible  to  determine  the  probability 


ul  e':i;  Itjtioil  dtlUVt-  <j  jtlO  il  -U  till:'  f'liiti  jbl  I  I  f  g  ul  f  t  iti  I  t  I  I.rll 

Into  the  continuum,  which  ic  incluijp.j  m  the  formei.  (The 
particularities  oi  the  distritiution  in  the  continuum  will  be 
investigated  in  another  paper). 

Several  characteristics  of  the  excitation  were  computed  by  the 
described  numerical  method.  Among  them,  the  most  important  were  the 
distribution  over  unperturbed  levels  moment 

M,  -  (<n>-no)/no  ,  the  2nd  moment  -  <i;Ari)">  =  <(n-<ri>)'^>.'  rig^.  and 
the  probability  of  excitation  to  high  levels.  In  order  to  describe  the 
latter  we  considered  the  probability  W,  ^  of  excitation  to  states  with 

n>  [l.Shoi  where  [  |  meahs  the  integral  part.  For  computations  in 
Sturm  base,  this  probability  included  also  the  probability  of  ionisation, 
namely  W'|  5  is  the  total  probability  in  states  r\±  1.5  no  plus  the 
probability  in  the  continuous  part  of  the  spectrum.  In  order  to 
eliminate  fluctuations,  we  also  determined  the  distribution  f^^  averaged 

over  At  periods  of  the  field;  as  a  rule,  At  was  chosen  40  or  60. 
Finally,  we  determined  the  average  distance  of  the  electron  from  the 
nucleus,  <2>. 

The  accuracy  of  the  numerical  results  was  checked  as  follows. 
First,  in  order  to  check  that  continuous  spectrum  was  being  properly 
taken  into  account,  we  performed  a  senes  of  experiments  with 
frequencies  larger  than  the  l -photon  ionization  threshold,  Wo  >  no/2. 
In  the  absence  of  resonances  within  the  discrete  spectrum,  the 
probability  on  discrete  levels  with  n  >  1.5  Oq  was  then  negligibly 
small,  so  that  the  probability  of  ionization  W,  ^  W,  5.  An  example  of 
dependence  of  W,  on  time  is  shown  in  Fig.  1.  in  Fig.  2  we  show  a 
comparison  of  the  theoretical  ionization  rate  with  the  numerically 
obtained  one.  As  can  be  seen,  there  is  an  excellent  agreement  with 
the  theory  of  l-photon  ionization  (2/),  which  indicates  that 
computations  in  Sturm  base  efficiently  rep'Oduce  continuum  effects. 

A  different  type  of  check  was  gotten  bu  increasing  the  number  L  of 
integration  steps  per  period.  The  relative  changes  of  ttie 
characteristics  of  excitation  produced  in  this  way  were  very  small. 
For  instance.  In  the  1st  method  (unperturbed  base.  iJEi)  with  n,-,  ^  66, 
<jOo  =  1-5,  to  =  0.04.  a  change  of  L  from  200  to  500  for  r  =  120  led  to 
a  relative  change  ^/w,  ^  =5  a<Z'-/<z>  -  5  I'v-*.  nr  the 
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irnjil  vjii.je':  ot  V/j  t.  For'  mstJnce.  tor  iIq  =  66,  -  2,  ~  O.O:.,  r  - 
M'lj,  Wj  c  -  '4. 1''"*.  Upon  rhjnging  L  from  100  to  200  the  relative 
change  in  proPatulitg  and  in  <2>  were  A<2>/<2>-:AW|  c/W|  c,  5. 10'^ 
We  can  therefore  assume  that  for  sufficiently  large  L  the  effects  of 
numerical  discreti2atton  in  the  integration  of  (30)  and  (34)  become 
negligibly  small,  and  have  no  influence  on  the  physics  of  the  problem. 

A  further  check  consisted  in  changing  the  total  number  of  levels 
both  in  the  bturrn  and  in  the  unperturbed  base,  and  also  in  matching 
the  excitation  characteristics  obtained  by  the  two  different  methods. 

One  such  comparison  is  shown  In  Fig.  3a.  where  it  can  he  seen  that 
there  is  a  good  agreement  between  results  of  computations  in  LiB  and 
in  5B,  and  also  that  an  increase  in  the  number  of  Sturm  levels  does  not 
change  significantly  the  excitation  probability  (which  Includes 
continuum). 

Such  an  agreement  not  only  takes  place  for  integrated  characteristics, 
but  also  for  the  disthbution  over  unperturbed  levels  (Fig.  4).  It  is  then 
possible  to  conclude  that  continuum  effects  do  not  lead  to  substantial 
modifications  of  the  excitation  dynamics,  at  least  for  not  too  strong 
fields  and  high  frequencies,  noreover,  the  Sturm  base  used  in  our 
computations  appears  large  enough  to  provide  a  satisfactory  model  for 
quantum  dynamics,  including  continuum. 

Of  course,  a  numerical  scheme  whatsoever  necessarily  involves  a 
discreti2ation  of  the  continuum,  and  shall  therefore  fail  under 
sufficiently  fine  tests.  Our  oven  method,  as  discussed  above,  correctly 
describes  the  continuum  spectrum  at  least  m  so  lar  as  one -phot  on 
effects  are  involved.  A  more  delicate  task  vv'ould  be.  lor  instance, 
reproducing  tunneling  in  a  static  field;  this  is  an  important  problem 
for  the  computer  simulation  of  actual  experiments  on  microwave 
lonication  as  we  shall  discuss  in  sec.  4.  Hern,  we  may  need  more 
sophisticated  techniques.  While  this  is  a  real  problem  for  future 
invest igai ions,  the  really  irnpoi  taut  question  nov:  iv  whether  our 
scheme  was  good  enough,  that  the  iocalication-delocalisation 
rnec tiani srn.  which  is  the  central  object  of  the  present  vvorl.  can  be 
considered  an  eiiective  phenomenon  and  not  jiist  an  artel  act  ot 
numerical  simulation.  In  this  respect,  the  agreerne.at  we  found  between 
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i.l jiiic-Jl  arid  qudiituiTi  cornfiutatiuns  in  the  delocalized  reqirne,  as.  we 
■ihall  discuss  in  the  next  section,  provides,  in  our  opinion,  the  most 
convincing  element  in  support  of  our  methods. 

The  dynamics  of  quantum  excitation  was  investigated  for  no  -  30, 
45,  66,  100,  and  the  field  ranged  in  the  interval  0.01  <  er,  <  0.34.  In 
order  to  facilitate  conversion  to  physical  units,  we  note  that  for  Oq  = 
100  the  frequency  =  lOGHz  corresponds  to  u>o  =  cuno^  =  1.51998, 
and  tn  -  ^rio'*  =  0.1  corresponds  to  €  =  5.14465  V/crn. 

For  clarity's  sake  we  have  grouped  our  numerical  results  following 
the  order  of  the  previous  theoretical  analysis.  Therefore  we  shall  now 
discuss,  in  turn,  the  results  on  the  classical  model,  the  results 
demonstrating  the  localization  phenomenon,  and  the  results 
illustrating  the  dependence  of  the  excitation  probability  on  the  field 
frequency. 
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4.2  fiumencal  results  on  the  Classical  Model. 

The  dependence  of  the  excitation  proPability  of  the  classical  system 
on  the  frequency  w,;,  and  intensity  Cg  of  the  field  is  shown  in  Fig.  5.  Here 
the  excitation  probability  W,  p,  is  computed  after  c=40'j^g  periods  of  the 

external  field,  v/e  recall  that  the  initial  value  is  irrelevant  due  to 
the  scaling  property  of  the  classical  rnotion.The  characteristic 
oscillations  with  minima  near  integer  values  of  Wy  are  connected  with 
the  presence  of  nonlinear  resonances,  the  strongest  of  which  correspond 
to  integer  ■■.Oii.  In  fact,  the  destruction  of  the  centers  of  resonance 
regions  occurs  for  larger  field  than  their  overlapping  (*). 

Then,  for  not  too  strong  fields,  a  part  of  the  trajectories  from  the 
initial  distribution.  vYhich  is  uniform  in  space,  upon  entering  the  chaotic 
region,  diffuse  to  higher  values  of  n,  but  the  rest  fall  into  the  central 
stable  region  of  resonance,  where  they  remain  giving  no  contribution  to 

The  characteristic  dip  for  cn'o  =  0.5,  which  was  also  observed  in 
numerical  experiments  on  2-dirn.  atoms  /'40/  corresponds  to  a  2nd  order 
(half-integer)  resonance.  The  sharp  maximum  of  W,  p  for  ■.ny  0.? 

(weakly  depending  on  Cy)  due  to  the  fact  that  for  this  frequency  most 
trajectories  fall  into  the  stochastic  layer  of  the  separatrix  of  the  big 
fundamental  resonance  Wy  =  1.  Already  after  half  a  turn  around  the 
resonance  they  pass  into  the  high-n  region,  where  excitation  is 
significantly  stronger.  An  analogous  excitation  mechanism,  connected 
Yvith  the  2nd  order  resonance  at  coy  =  0.5.  explains  also  the  maximum  at 
(Oy  =  0.43. 

In  the  classical  system  diffusive  excitation  takes  place  only  when  the 
field  strength  exceeds  the  critical  value  for  which  the  last  KAM 


(*)The  value':,  of  the  field  at  w'hich  tlie  centers  of  resun.ance  reqions  are 
de-ftioyed  •were  numencally  deteirmiied  in  ."I"?,-'''. 
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inv.in.jnt  iurve  li  iJestrciyniJ  dud  thert-  it.  a  frjnt.if lOfi  to  qIoD'jI 
■tUocasticity  (see  e.g..  /5/'i.  From  Fig.  5  we  see  that  the  actual  value  of 
tor  (ji>o'^  1  IS  near  to  0.02,  which  satisfactorily  agrees  with  the 
theoretical  value  (3)  obtained  by  the  resonance  overlap  criterion  /5/. 

Fig.  5  gives  an  overall  idea  of  the  classical  behaviour.  Other  numerical 
results  such  as  the  comparison  with  the  solution  of  the  diffusion 
equation  or  with  the  quantum  distribution  on  the  unperturbed  levels  will 
be  given  in  the  following  sections. 

4.3  The  Distribution  over  the  Unperturbed  Levels 

Here  we  shall  describe  the  features  of  the  numerically  computed 
quantum  distribution  over  the  unperturbed  levels  in  the  various 
parameter  regions  which  have  been  discussed  in  our  previous  theoretical 
analysis.  In  this  way  we  shall  show  that  numerical  results  support  the 
theoretical  estimates  given  above. 

For  high  levels  (  e.q.  ho  ^  100  )  and  tQ  -s  ^  0.02  the  perturbation 
strength  V  =  (3/2)  n^c  is  significantly  larger  than  the  level  separation: 
V/AE  ~  (3/2)  €o  h  ^  •>  so  that  the  field  would  be  expected  to  connect  a 
number  of  unperturbed  levels.  Yet,  even  for  Cnh  >>  l  no  diffusive 
excitation  will  be  observed  if  €<(^^  .  (The  opposite  case  Cci  n 

corresponds  to  the  region  below  the  "quantum  stability  border"  /35/). 
This  is  illustrated  in  Fig.  6,  where  an  example  of  stationary  distribution 
in  the  region  of  stability  is  shown.  This  distribution  remains 
essentially  unchanged  upon  further  increasing  the  computation  time. 
Classically,  this  fact  is  due  to  the  stability  of  the  motion,  and  qu-antum 
mechanically  to  the  very  small  probability  of  tunneling  into  regions 
classically  forbidden  by  smooth  invariant  curve  (oec.  2.?).  Howe/er,  for 
a  reliable  detection  of  the  tunneling  descnoed  in  Sec.  2.5  particularly 
accurate  investigations  are  required.  It  is  also  desirafiie  to  increase  n,)  , 
because  even  for  ry,-  lOO  the  tunneling  region  appears  rather  narrow  (see 
eq.27).  Nevertheless,  we  thinF  that  tunneling  excitation  can  be 
investigated  both  in  numerical  and  in  laboratory  experiments,  where  at 
the  present  time  it  is  possible  to  prepare  states  Vvith  iiq  /25/.  We 
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field'::  are  e-till  lari  irnj. 

For  Held  srrength  exceeding  tne  cnncal  'calue  <7).  difiui.r.'e  excironrui 
takes  place  in  the  classical  system.  However,  in  the  quantum  case,  for 
field  strength  lesser  than  the  delocalisation  border  ilS)  the  phenomenon 
of  quantum  localisation  is  observed,  in  consequence  of  which  the 
distribution  over  the  unperturbed  levels  reaches  the  stationary  lorrn 
(lo)  and  then  does  not  change  upon  increasing  the  time  of  interaction 
with  the  field.  In  this  situation,  the  ionization  probability  is  very  small, 
and  can  be  neglected  tor  tne  given  interaction  time.  A  typical  example  oi 
guanrurn  localization  is  shown  in  Fig.  ?.  Here  we  see  that  classically 
there  is  a  diffusive  excitation,  so  that  the  classical  distribution 
obtained  by  the  numerical  simulation  'satisfactorily  agrees  with  the 
theoretical  formula  (5a)  (the  classical  border  of  stability  was  here 
chosen  at  FT  -  55  according  to  numerical  results).  The  quantum 
distribution  was  obtained  by  tlie  -riturm  base  rnettiod  with  H5  z  576;  here, 
as  well  as  in  Fig.  d,  there  is  a  good  agreement  with  the  results  of 
computations  by  the  unperturbed  base  method.  In  contrast  with  the 
classical  result,  in  the  quantum  case  an  exponential  drop  lollowed  by  a 
rnultiphoton  plateau  is  observed,  almost  unchanged  under  a  change  of  t 
from  120  to  600.  The  quantum  limitation  of  chaos  also  led  to  a 
significantly  lesser  excitation  probability  m  the  quantum  than  in  ttie 
classical  case  (see  also  Fig.  3). 

Another  convincing  manifestation  of  localization  was  the  saturation  of 
the  diffusive  growth  oi  the  moment  oi  the  quantum  distribution  (Fig.  6). 
The  agreement  between  quantum  and  classical  dynamics  here  holds  only 
over  a  small  initial  time  interval  The  smallness  oi  fp,  is  due  to 

the  smallness  oi  the  cias-.iral  dniusion  rate, 

I'l  Fig.  'J  It  IS  hewn  hu;-.  the  n"iiiiolised  j  ei  aye  CwtauLf 
RL - <z( t )>.  11,1^  of  the  electron  from  the  nucleus  depends  on  tinie.  In  the 
i.l’i,.s)i  ,-ji  .ypteiFi  Ihis  ill  .t-ifii  »'  i|tiiw  irid  thi^  hIhi  linii  I.ar  iv.  ii.j 

irorn  tne  nucleus.  Insfead.  m  the  yuantijr''  <:ase.  nwinq  tn  gu-jntum 
I  O'  al  I  .'at  1 1  III,  the  e|ei  tti'iii  1  eep  ii  s  1 1 1  at  Uiy  aicmnd  Its  initial  pnsitii;ih. 

I  he  iei)irne  lu  quahtum  h.n. ali.~ai loti  v-.a'.  ihvestiqated  im  initial  levels 
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ii,,  -  I,  45,  66, 1 UO,  lot  I i‘eC|ut'iiC  It-:-  I  ■  -'3  and  Held  interiiltleS  U.03 

to  ^  0.12,  and  localization  length  l>>l  (which  is  the  condition  for 
applicability  of  the  estimate  (ID).  Numerical  values  for  the  localization 
length  were  determined  directly  from  the  stationary  quantum 
distribution.  Comparison  of  these  numerical  data  with  the  theoretical 
values  (  Eqs.  13,14)  yields  good  agreement  (Fig.  10).  The  observed 
dispersion  of  points  is  apparently  connected  with  the  presence  of  islands 
of  stability  in  the  classical  system.  Also,  the  presence  of  a  typically 
quantum  resonance  structure  niay  play  a  role  in  this  respect. 


If  the  field  exceeds  the  quantum  delocalization  border  (15)  then  the 
quasi  classical  diffusion  over  the  levels  is  suiiiriently  fast  and  no 
localization  takes  place.  In  this  regime  the  evolution  of  the  distribution 
function  can  be  approximately  described  by  the  diffusion  equation  (4).  An 
example  of  distribution  in  the  delocalization  regiirie  is  shown  in  Fig. 
1  la,  where  it  can  be  seen  that  the  quantum  distribution  agrees  with  the 
solution  (5a)  of  the  Fokker-Planck  equation. 

In  order  to  check  the  valididy  of  the  estimate  (15)  for  the 
delocalization  border,  we  investigated  the  dependence  of  the  excitation 
probability  on  for  different  values  of  rio  and  of  a'f,.  This  dependence 
on  the  rescaled  field  c®  =  shown  in  Fig.  13.  For  each  value  of 

no,  (Oo,  the  excitation  probability  W,  ^  was  also  rescaled  to  the 
corresponding  classical  value  taken  for  Cq  -  'q*''’-  other  words,  in  Fig. 
13  we  actually  plotted  W|_5(f~o)  =  =  £q''*- 

In  case  that  delocalization  should  actually  take  place  for  and 

that  in  the  delocalization  regime  the  excitation  probability  should  keep 
close  to  Its  classical  value,  then  all  the  lines  showing  the  dependence  oi 
W]  c  (tn)  (or  different  values  of  n,i,'j.',-|  would  be  expected  tn  meet  tor  Ci-,  = 

1  at  the  value  Vi,  c  -  l.  As  can  be  seen  from  Fig.  13.  this  is  pjst  what 
actually  happens. 

An  interesting  feature  of  Fig.  13  is  that  the  dependence  of  the 
ionization  probability  on  the  field  strength  at  fixed  n,;,.  0.1^  is  not  always 
monotonic.  For  example,  the  data  corresponding  to  n,i  =66,  mq=3  clearly 
indicate  a  'bump"  occurring  in  the  ionization  curve.  The  existence  of 
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lecently  pomteJ  out  in  .-M?/  -and  j  thenrHtir.al  e/planation  wa--:-  put  forth 
tn  /4i:i/. 

In  the  localization  regime,  the  dependence  of  the  excitat'on  proPaPilty 
on  field  intensity  can  be  approximately  described  by  V/,  c  cc  tf/K  Fig. 
13  also  clearly  indicates  that  the  experimental  value  of  f  changes 
substantially  with  n,,,  cuf,.  so  that  loining  oi  all  lines  at  a  single  point 
for  tn  =  1  IS  nut  a  trivial  occurrence,  and  can  be  considered  as  a 
confirmation  for  our  estimate  (15).  This  diversity  in  the  values  of  f  is 
connected  with  the  different  number  of  photons  which  are  required  for 
excitation  in  states  with  n  >  [l.5  n,;,].  However,  the  experimentally 
determined  value  Kg  is.  typically,  substantially  lesser  than  the  number  Kg, 

of  photons  theoretically  reguired  for  direct  transition  from  rin  to  n  1.5 
rifl.  which  is  kg,  =  I  ♦  [Siin/IHcoo).  For  instance,  for  nr,  -  66.  coq  =  1.  one 
has  i  g  ~  7.  kg,  =  19. 

In  our  opinion  this  difference  is  due  to  two  effects.  The  first  is  that 
rnultiphoton  transitions  do  not  necessarily  start  from  the  initial 
unperturbed  state,  but  may  start  from  anywhere  inside  the  stationary 
distribution  (10)  which  sets  up  after  a  while.  In  other  words,  when 
1  •••  1  excitation  may  start  from  levels  n~  nrj+  1.  and  this  reduces  the 

multiphotomc  degree  k.  The  other  reason  is  the  appearance,  for  high 
levels  n>ni-i,  of  a  rnultiphoton  plateau  of  equidistant  resonances  /2I/. 
Examples  of  distributions  7^  which  clearly  exhibit  this  rnultiphoton 
plateau  are  given  in  Fig.  14a,b  (see,  also  Fig.  4).  The  differences  in 
unperturbed  energies  =  -  l/2n‘  between  consecutive  peaks  of  the 

distribution  are  equal  to  the  field  frequency:  therefore,  the  sequence  of 
pejls  can  be  naturally  explained  as  the  result  cf  a  chain  of  one-nhoton 
transitions. 

In  ftie  I  ases  illustrated  by  Figs.  Ha,  b  thes,;'  transitions  start  directly 
from  the  initial  state  ng,.  and  the  peal  s  can  be  enumerated  simply  by  the 
number  oi  photuns.  Howevc: .  the  situalion  is  not  always  that  simple;  in 
1  ceiies  oi  cases,  the  chain  oi  peals  does  not  start  iroin  rif,.  but  rather 
fiuni  sornevvheie  inside  the  localized  distribution  (see  Fig.  j  in  Ref. 
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/21/,',  oiiil  it  li  even  possible  to  observe  two  or  tlii'ee  distinct  clioins 
within  the  same  distribution. 

On  the  high  levels  the  amplitudes  of  the  peaks  become  roughly  the 
same  and  they  build  up  an  equidistant  plateau.  Increasing  n  still  further, 
the  peak  amplitudes  do  not  decrease;  this  seems  to  be  due  to  the  fact 
that  on  high  levels  the  field  is  strong  enough  for  the  probability  of 
transitions  between  nearby  peaks  to  be  significant  (saturated 
transitions).  This  is  the  second  reason  why 

Upon  increasing  the  field,  the  multiphoton  plateau  rises  as  a  whole 
(Fig.  H).  The  resonant  peaks  become  broader,  but  in  a  number  of  cases 
they  do  not  disappear,  even  in  very  strong  fields  and  in  the 
delocalieation  region.  However,  this  can  usually  take  place  only  for  large 
<Oo  (compare  delocalization  in  Fig.  |  la  and  Hb,  for  Oq  =  tOO, 
o.>o  =  1.5  and  3,  respectively). 


In  our  opinion,  the  appearance  of  the  multiphoton  plateau  below  the 
ionization  threshold  is  in  its  substance  akin  to  the  appearance  of  peaks 
in  the  energy  distribution  of  photoelectrons  which  is  observed  above  the 
ionization  threshold  /41/42/.  Indeed,  for  large  n  the  distance  between 
nearby  levels  is  very  small:  and  the  spectrum  in  this  region 
behaves  like  a  quasi-continuum.  It  is  then  reasonable  to  expect  that  the 
peak  structure  observed  in  the  discrete  part  of  the  spectrum  will  persist 
also  in  the  continuum. 

It  is  possible  that  a  theoretical  explanation  of  the  multiphotonic 
plateau  in  the  under  threshold  distribution  may  be  given,  along  similar 
lines  as  in  /43/. 
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All  f.-.jiiiplf  iji  deptriderice  Of  'w'j  c  on  ‘jL'n  for  fixed  c,-,  =  u.OA  ond 

different  iiu  is  shown  in  Fig.  15.  It  is  hene  opporent  thot  for  high 
It s^gijencie:-.  (o,,  >  i  the  e.':cit.jtiori  proboDintg  is  signinooritlg  less 
I’ll  111  Mi*^  sorresponding  ciossic  ji  one.  I  he  reoson  is  Uiot  ffie 
deloc ol  1  sotion  threshold  '15'  increoses  with  the  freguencg  ovi.  so  that 
3  irijpjritg  of  points  in  Fig.  15  belongs  to  the  region  of  localisation. 
Vv'e  recall,  however,  that  the  estimate  (15/  is  onlg  valid  lor  con  :■  i. 
Indeed,  tor  ■..■a  I  a  dgnarnical  amplification  or  the  classical 
e'v.itation  tales  place,  as  described  in  section  7. A  isee  Fio.  5),  From 
Fig-.  5.  15  it  appears  that  for  rrj  -  c'.c'd  the  mavamuni  classical 
..i'lfation  IS  gotten  tor  cOi-,  -a  n.,'.  (-jci'A'.  as  is  seen  trorn  tig.  i -.  tne 
guariti.irn  probabilitg  of  excitation  for  n,-,  -  45.  66.  icici  i;,  close  to 

'’he  classical  value,  and  even  the  corresponding  distributions  on  levels 
'Fig.  16/  loof  rather  close  to  classical  results,  'w'e  niterpret  this  last 
excitation  in  the  sense  that  in  the  guantum  sgstern  lor  thi-  ireguencg 
there  is  a  delocalization.  Thus  numerical  experiment:  :tiow  that  I'jc 
sufficientig  weal  field  <■.,)  -  'j.';'-i  delocalization  c 'n  ‘aie  place  at  a 
Ireguencg  coi-i  t-  i.i,,',  lovver  tnan  tne  t  epier  treguenc'j. 

Ttie  line  structure  ot  the  dependence  oi  excitation  on  ireguencg  ic 
shown  in  Fig.  17.  In  the  localization  region  ■■.■•i-i  i  nnr.  ripserves  an 
essentiallg  resonant  dependence  on  freguencg.  Fm  low  treguencies 

•=■  'o.c  most  resonances  disappear  and  the  dependence  on  '■j,',',  becomes 
■■■■moother.  An  analogous  smoothing  occurs  in  tne  region  -t  i,  upon 
increasing  rin  from  70  to  100. 

For  still  lower  ireguencg,  'V',i  c'.r',  i/ne  i.ill,  '.nfr:  th.'  o 

classical  stabilitg:  therefore,  excitation  i.x/  rtw  ■:l'i  /.i  jl 

jgstern.  Then,  tlie  guanf/im  excitation  :naiph_i  iiiir.itii  !:>■:,  Im.i.  Inn  , 
the  dependence  of  the  excitation  probabilitg  cn  trequenc/j  has  a 
thi  estiold  rtiaracter,  in  that  lonication  t.jif;.  niaro  'oi 
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Till'..  I."';.!  iiti.jie  for  the  i;h.30tjc  threT-hold  <x>^  io  luotifipij  bu  the  f.3i:t 
that  tor  'Oy  I  there  are  no  I'jt  order  resonances  between  the 
frequency  ot  the  external  field  and  the  harmonics  of  the  frequencies  of 
the  motion  of  the  classical  electron  /1 7/;  therefore,  chaotic 
excitation  for  frequency  a'o  <  I  can  take  place  only  for  so  strong 
tields  that  hiqher-order  resonances  overlap.  On  the  other  hand,  for 
very  low  frequency  cOf,  =<  l  the  value  of  the  critical  field  c,  coincides 
with  that  for  the  classical  static  field  ionization.  Notice  that  the 
threshold  (36)  holds  only  if  Cq  >  1/50.  Otherwise,  the  chaotic 
threshold  has  to  be  determined  from  eq.  (3),  and  is  equal  to  (jl>/:I/ 

(50(0^. 

According  to  the  theoretical  estimate  (2*1).  the  diffusive  ionization 
IS  more  eifective  than  direct  one-photon  ionization.  In  order  to  check 
fhis  prediction  we  performed  a  series  ot  numerical  experiments,  in 
which  the  ionization  probability  from  states  with  riu  =  30  or  rip  =  66 
was  investigated  over  a  broad  range  of  frequencies.  The  field  intensity 
was  so  chosen,  that  direct  two-photon  ionization  was  considerably 
lesser  than  I -photon  ionization  except  for  intermediate  resonances; 
moreover.  Cn  <  Cj,  where  e.  is  the  critical  intensity  for  static  field 
ionization.  In  this  situation,  the  photoeffect  is  expected  to  display  a 
threshold  dependence  on  frequency,  with  negligibly  small  ionization 
probability  for  •j.Vj  =  nci/2. 

5uch  a  picture  of  the  photoeffect  proved  to  be  incorrect.  In  Figs.  16, 
19  we  show  the  dependence  on  frequency  of  the  excitation  probabilty 

above  a  level  FT  after  a  dimensionless  time  t  =  Tcu.On  which 
corresponds  to  the  same  physical  time  for  all  'ji'n’s.  Computations  were 
made  in  Sturm  base,  so  that  encloses  the  probability  of  iransition 
into  the  continuous  spectrum.  For  n,;,  =  FT  -  9'):  therefore  m  thr:, 
case  Wj  ^  W|.  For  On  =  66,  we  took  iT  =  99,  and  then  enclosed  a 

significant  part  of  the  probability  on  discrete  unperturbed  le"els.  In 
Figs.  Id,  19  the  most  effective  excitation  is  observed  at  irequencies 
'well  below  the  i  -  fihoton  threshold.  The  ne'w  ihre'  iiold  vilue  o.i .  •  i  is 

close  to  the  correspondinq  classical  value  and.  as  eiiplamed  above,  is 


■  JrfHf  lUlflPil  tiq  Utf'  C'.Kidl  i  lijll  o!  u'-'i:'!  i  J(i;i  1 1 IIJ  III  .'I  IJ  '.I' |J*'I  I  I?  :  i5M  il  |i;  P  ;  . 

Flu  :uf iiL'ietitlq  L'ti'onq  lielJ.  ■.o,  jppeoi';-  L'lgnitu: jhUm  leLTL'i  ttijci  I: 
inr  'll  :  L'.  1.  n,-,  ;  T'O,  the  PiLpenment.jl  vilue  for  ''jv  ip  '■  inq  for 

i;i'i  =  u.Oci,  ri,-|  -  66  the  thresholLi  is  u.'i5. 

In  the  QUjhn.irn  ojse.  however,  there  is  stiii  .3  snijll  proDoDility  01 
e 'LL itdtiori  even  from  the  region  'jl'i,-,  <  0'!.,  diie  to  tunneling  into  the 

L'hj'S'sicjlly  forbidden  region.  In  the  interval  'x'^  <  ix'ij  '.1.^  = 
delocdliratiDn  toPe'S  place,  so  that  is  close  to  its  classical  value. 
For  i'.0|_  r  'On  ''O4  guanturn  ettects  lead  to  the  localization  ot 

diffusion:  therefore,  in  that  region  one  observes  an  essentiallg 
different  ericitation  picture  than  in  the  class-ical  system.  A  decrease 
01  VvV  'With  (.Oij  IS  observed  also  in  the  classical  case,  due  to  the 
decrease  of  the  diffusion  rate  'vvith  ••on  (see  eg.  (5)1.  Indeed,  even  in 
the  presence  of  diffusion,  the  excitation  probability  above  a  given 
level  n  is  nearly  zero  if  at  the  given  observation  time  tg  the  diffusion 
has  not  yet  reached  this  n  (i.e.  if  <<  in-nn)‘-/[c'.  Ho'wever.  in  the 

guanturn  case,  isolated  ionization  'Spifes  survive  in  this  treguency 
region  due  to  transitions  into  the  continuum  ■■■•la  interinediate  resonant 
levels. 

In  the  freL|i.iency  region  otri  •.i.Vj.  tlie  probability  sharply 
increases,  because  01  the  possibility  of  direct  one-photon  transitions. 
Ho'vvever,  in  this  region  is  significantly  lesser  than  tor  ■'.I'lij  •'I'l. 
■when  diffusive  ionization  occurs.  In  the  region  I'x'ii  >  '.ilj,  nuinerical  data 

■satisfactorily  agree  with  the  theoretical  formula  il'J).  Tins  agreement 
indicates  that  the  .riturrn  base  rnethod  01  integration  eiiiciently 
iiescriDe;  continuous  spectrurn  eiiect':-.  6ee  also  tig.  ',190). 

in  fig.  19  vve  see  tbat  loi  ‘■.l-■'||  ".'is  the  gnantum  protiability  lU 

e.’icitation  is  close  to  the  classical  one;  again,  as  in  the  case  of  Fig. 
1  j,  fhi’  (oean:.  that  in  this  ireguency  region,  for  n,)  -  06,  tr,  =  "."d, 
'her''  !.  guanturn  delocalization. 

'  .  Ini' 


iinini'rii  ai  I'Zpi'rimeiit';.  indicafe  tnaf  the  tnresnolo  v  ilue 


'JnfuMVi?  I tJtiun  in  thi:-  trHqut’ficy  reqion  lii?:.  bc'twt'f'ri  v'.'.M  jnd 
'j.05.  (Fig.  5).  In  the  quantum  case,  upon  vanjjnq  Cg  from  0.0-4  to  (uOS 
the  excitation  prohatnlity  changes  by  about  two  orders  of  magnitude. 
This  means  that  delocalization  occurs  already  at  -3  so  that  the 

field  intensity  yielding  10.^  ionization  (  which  was  studied  in  ref.  (26) 
will  be  found  to  agree  with  classical  predictions.  In  other  words, 
these  numerical  experiments  ot  ours  snow  that  laboratory  experiments 
/ 12/26.''  were  performed  in  that  parameter  region  where  l-dirn. 
delocalization  (and  a  fortiori  2  dim.  delocalization)  takes  place,  and 
this  explains  the  observed  agreement  with  predictions  from  the 
classical  model. 


■4.5  Stability  of  Quantum  Oiffusion 


Even  though  the  "diffusive"  ionization,  taking  place  in  the 
delocalized  regime,  is  to  some  extent  similar  to  the  classical 
diffusion  which  occurs  in  the  chaotic  regime,  the  quantum  system  is 
still  short  of  exhibiting  all  the  statistical  properties  that  would  be 
expected  of  classical  chaos. 

The  most  striking  difference  is  the  absence,  in  quantum  dynamics, 
of  the  strong  instability  and  of  the  rapid  loss  of  memory  associated 
with  classical  chaos.  In  computer  experiments  this  effect  leads  to 
irreversibility.  Indeed,  even  though  the  exact  equations  of  motion  are 
reversible,  nevertheless  any,  however  small,  imprecision  in  solving 
them,  such  as.  e.g..  computer  round-off  errors,  is  magnified  by 
exponential  instability  of  orbits  to  the  extent  that  initial  conditions 
are  eitaced  and  reversibility  is  'I'lereiore  destroyed. 

Investigations  aimed  at  veriiyinq  whetlier  an  analuyous 
irreversibility  would  be  displayed  also  by  the  numerically  computed 
i;|uantnrn  evolutnifi  vvere  described  m  /6/  tor  the  mcled  roiaior.  Here 
w'e  V'l'ill  present  nutneric.al  resnits  for  time-reversal  e^'periments  on 
ttie  1  -  dimensional  H-afom.  (Fiys.  II.  12'.  Th>’  chiisen  p irairoter 
'■'alues  lie  in  the  region  of  deloc .aliz.atinn;  therefore,  up  to  ttie  mornent 
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■jf  t iin>?- ir  -  6'?')  iliin.rE.ive  1 1 of toif  if.  qouiQ  on,  Doth  in  the 

qu.jntum  .3niJ  in  ttie  classic.sl  -system.  Indeed,  the  di-stribution  on  the 
quonturn  levels  .at  r;  =  dO  is  close  to  the  cl.assical  one,  .and  is  ■■.■veil 
descntied  by  tormul.a  t  5a).  Then,  at  t  =  60  -//e  reversed  the  velocities 
of  all  particles  (N  -  1000)  in  the  classical  ensemble,  and  chanqed  the 
■*-ave  fiihction  of  the  quantum  atom  to  its  comple:-;  coniuyate.  In  both 
classical  and  quantum  mechanics,  the  H-atorn  ■v'could  be  expected  to 
find  its  'vvay  back  to  the  initial  state.  Ho-Vv-ever,  due  to  the  finite 
computer  precision,  in  the  classical  case  such  a  return  is  not 
observed.  The  system  retraces  back-.Yard  its  history  just  for  a  fe-'V 
periods  of  the  field,  and  then,  again,  diffusive  excitation  occurs. 

Instead,  in  the  quantum  case  an  almost  exact  reversion  of  motion  is 
gotten;  at  time  t  -  120  the  electron  comes  back  to  the  initial  level. 
This  IS  even  more  remarkable  on  account  of  the  fact  that,  in  order  to 
restore  the  initial  state,  some  of  the  total  probability  had  to  be  called 
back  from  the  continuum. 

The  conclusion  must  be  drawn  from  this  exact  reversibility,  that 
even  though  Che  quantum  diffusion  wivcli  occurs  in  the  delocalized 
regime  of  the  H-atorn  is  by  no-vv  the  most  chaotic  example  of  quantum 
motion  hitherto  investigated,  nevertheless  this  quantum  "chaos"  is 
essentially  different  from  the  real  chaos  of  classical  dynamics. 


k 


5.  Experimental  Results 


A  large  numDer  ot  (aooratory  expen  men  rs  on  nyoroqen  ano  air  an 
atoms  in  highly  excited  states  have  been  performed  up  to  now 
/1 2/ 1 3/22-26/.  Additional  interest  for  such  experiments  has  recently 
arisen  m  connection  with  the  possibility  of  chaotic  motion  in  quantum 
mechanics.  It  is  now  possible  to  perform  experiments  on  microwave 
ionization  on  atoms  prepared  in  extended  quasi  I -dim.  states  /23/.  Here 
it  is  possible  to  measure  the  ionization  probability  identified  with  the 
population  of  levels  higher  than  some  suificiently  large  n,  including 
rontinuum:  this  definition  is  particularly  convenient  lor  comparison  with 
numerical  experiments.  Also,  it  Is  possible  to  measure  the  probability 
distribution  on  unperturbed  levels.  This  allows,  in  pnncipie,  for  a 
careful  comparison  of  experimental  and  numerical  data. 


In  particular,  very  accurate  experiments  on  I -dim.  H-atorns  were 
carried  through  as  described  in  /23/.  The  range  of  parameters  for  these 
experiments  lies  inside  the  region  of  low  frequency  (Wf,  s:  0.2)  and  of 
classical  stability,  so  that  the  results  can  not  be  used  as  a  test  lor  the 
theory  presented  in  this  paper. 

A  different  series  of  experiments  /26/  was  performed  on  2-d 
H-atoms.  The  conditions  of  these  experiments  not  only  lie  above  the 
2-dim.  delocalization  border,  but  even  above  the  l-dirn.  one.  For  this 
reason,  our  results  predict  an  agreement  with  classical  computations,  as 
indeed  was  found  in  /26/.  One  possible  explanation  for  the  not  complete 
agreement  obtained  is  that  the  experimental  values  in  /Zb/  migtn  be 
above  but  close  to  the  delocalization  border,  when  one  should  not  expect 
a  better  agreement  than  within  a  factor  2.  To  clarify  this  point  we  show 
in  Fig.  zu  the  comparison  ot  ntirrieiir ai  i-di,n.  quantum  .and  classical 
ionization  probability  slightly  above  the  border:  there  is  a  strong 
excitation  in  both  cases,  but,  unliFe  strongly  delocalized  cases.  (Fig. l  la! 
here  the  two  results  only  agree  '.within  a  Suiti. 


jl  tSLhMiQut'S  jllow  lor  vt'ru  jco.irjto  ii'ic'jrurornt'tit.;., 
it  I'i;  highly  desirable  that  the  conditions  of  the  experiments  be  defined 
as  precisely  as  possible;  for  example  it  is  more  convenient  to  choose  a 
single  excited  state  than  a  rtucrocanonicai  distribution,  from  our  point 
of  'uew,  a  most  important  goal  for  future  experiments  is  to  observe  and 
to  study  the  new  and  unexpected  localization  phenomenon  in  classically 
chaotic  situations.  For  this  it  is  necessary  that  the  frequency  (o,,  be 
increased  above  a'd  s  l,  since  in  the  region  a  large  separation 

between  the  classical  chaotic  threshold  'jli^  and  the  quantum 
delocalization  border  (o,  is  expected.  In  the  high  frequency  region  it  is 

also  possible,  by  varying  the  field  strength,  to  observe  the  transition  to 
delocalization  as  well  as  the  other  phenomena  aescribeo  in  the  present 
paper. 

Also,  m  order  to  give  experimental  evidence  for  the  "freezing"  of  the 
rvave  paci'et  in  localization,  it  would  be  desirable  to  dispose  of  a  control 
on  the  interaction  time.  This  latter  possibility  lies  within  the 
capabilities  of  present  day  technique.c'24/. 

we  would  lif  e  also  to  stress  that  all  the  phenorneha  described  in  this 
paper  should  be  observable  not  only  in  H-atoms,  but  also  in  different 
alFali  atoms.  In  order  to  produce  hydrogen-like  states  in  such  atoms,  one 
should  take  into  account  that,  the  unperturbed  spectrum  for 
highly-excited  alkali  atoms  is  slightly  different  than  in  H-atoms,  due  to 
quantum  defects.  However,  for  values  of  1>3  this  quantum  defect  is 

negligible.  Since  in  linearly  polarized  fields  the  magnetic  quantum 

number  rn  is  a  constant  oi  the  motion,  by  exciting  states  with  m  >  3  it 
IS  possible  to  excite  states  with  |>T>,  which  correspond  very  well  to  the 
hydrogeniL  situation.  It  is  then  possible  to  consider  also  1 -dimensional 

states  by  exciting  levels  vvith  m>T>,  ry.  =  o.  n|  -  n'|m|-l  and  this 

excitation  can  be  achieved  via  light-induced  resonant  transitions.  One 
would  then  get  a  situation  in  which  localization  and  other  effects  of 
guanium  chaos  might  be  studied. 


FIGURES  CAPTION 


Fig.  I  Ionization  probability  W,  as  a  function  of  time  z  (number  of 

microwave  periods)  for  the  case  n^  =30,  u)o=  30,  if,-  O.OTS.  The  solid 
line  is  drawn  accordtnq  to  the  analytical  expression  (22)  while  the 
crosses  are  the  results  of  our  quantum  numerical  computations.  The 
excellent  agreement  with  the  theory  even  for  very  large  frequencies  is  a 
check  01  our  numerical  computations  and  snows  that  Sturm  base 
efficiently  takes  into  account  the  continuous  spectrum. 

Fig.  2  Ionization  rate  versus  field  intensity  for  the  case  nu=30,  (i>o=30. 
Like  in  the  previous  fig.  1,  the  straight  line  is  drawn  according  to  the 
theoretical  expression  (22)  and  the  crosses  are  results  of  quantum 
numerical  computations.  Here  also  notice  the  very  good  agreement 
between  theory  and  numerical  results. 

Fig.  3a  Excitation  probability  W,  ^as  a  function  of  time  r  for  the  case 

nci=  DO.  to=0.0‘f,  <jOo=2.5.  The  quantum  numerical  computations  are 
performed  by  using:  i)  the  unperturbed  base  witn  N:I92  basis  set 
eigenstates  (  solid  line);  n)  the  Sturm  base  with  N=  364  (dotted  line); 
iii)  the  Sturm  base  with  N=  576  (dashed  line).  The  fairly  good  agreement 
of  the  three  curves  is  a  check  of  the  numerical  computations.  The 
classical  ionization  curve  is  also  shown. 

Fig.  3b  Classical  (1)  and  quantum  (2)  excitation  probability  W,  ^  as  a 
function  of  time  for  the  same  case  as  in  Fig.  3a. 

Fig.  4  Quantum  probability  distributions  Tin)  over  the  unperturbed  states 
averaged  over  60  values  of  c=ci.'t/2rf  within  tlie  interval  60'-'r<l20. 
Here  00=66.  i!:o=0.04.  'Co^l.S.  Three  different  curves  are  plotted 
corresponding  to  integration  in  Sturm  base  with  a  basis  set  of  N=  364 
and  N=576  basis  functions  and  integration  in  unperturbed  base  with  N= 
192.  The  three  curves  are  so  close  that  are  not  resolved  in  the  graph  and 
this  is  an  additional  check  on  the  accuracy  of  numerical  computations. 


Fi'^.  5  Llriiic  ol  lofa^atioi'i  tirobatiiliti^  Wj  c,  at  ter  c=-40  a»q  j;.  o  function 
of  the  microwave  frequency  for  different  microwave  intensities.  Here 
W,  5  IS  the  total  probability  above  the  action  value  n  ^^l-SriQ,  (x)  6g 
=0.02;  (.)  Co  (o)  (q  =  0.04;  (2i  )  -  0.05;  (■)  €(,  =0.06. 

Fig.  6  Classical  (-)  and  quantumt-i  probability  distribution  Hn)  averaged 
over  40  values  of  t  within  the  interval  60<r<120.  Here  05=  100, 
C(;i=0.01,  rofrl.S.  For  these  parameter  values,  cc<c^.|_-r^  and  therefore  both 
classical  and  quantum  packets  are  localized.  Notice  the  small  tunneling 
through  the  classical  KAM  invariant  curves. 

Fig.  7  Classical  (  dotted  curve)  and  quantum  (solid  curve)  probability 
distribution  Hn)  averaged  over  40  periods  of  r  for  the  case  riii  =  66,  aig  - 
2.5.  to  =0.04.  Fig.  7a  average  within  the  interval  60<t<120;  Fig.  ?b 
average  within  the  interval  5DO<t<DOO.  The  dashed  line  in  both  figures 
represent  the  analytical  solution  (5a)  of  the  Fokker-Planck  equation 
which  fairly  agree  with  the  classical  numerical  results.  On  the  contrary, 
the  quantum  distribution  is  localized  and  do  not  change  significantly  by 
increasing  the  interaction  time  with  the  microwave  from  120  to  600. 
The  only  difference,  as  expected,  is  the  slight  increases  in  the  peaks  of 
the  small  multiphoton  plateau. 

Fig.  8  Second  moment  hi:.  =<(n-<n>)‘>/no‘  of  the  classical  (solid  line) 
and  quantum  (dashed  line)  distribution  as  function  of  time  t=coT/2Ti  tor 
the  same  paremeters  of  fig,  7.  The  localization  of  the  quantum  packet 
sho'A'n  in  the  previous  fiq.7  leads  here  to  the  suppression  of  the 
diffusive  growth  of  the  rnornent  H_,. 

Fig.  9  Mormalized  average  distance  Rl  =  of  the  election  from 

the  nucleus  as  a  function  of  c  for  the  same  case  as  in  Fig.  7.  (dotted 
line)  r|u,3rir.um  case;  (full  line)  classical  case.  Also  here  the  quantum 
suppression  of  du fusion  is  clearly  manifest. 

Fig.  10  Localization  length  as  a  lunction  of  field  mtensity  for  diuerent 
parameters  value--.  The  dots  tfini-s.ppoiid  lu  niime:  icaliy  measured  values 


of  I  which  are  in  good  ogreernerit  with  the  solid  curve  given  hy  the 
analytical  estimate  ( 1 3) 

Fig.  1  la  Classical  (  )  and  quantum  (  )  distribution  function  Hn)  averaged 
over  40  values  of  t  in  the  interval  40  <r<30.  Here  hQ-  100,  a'n-l.S,  eo- 
0.08.  Notice  the  fairly  good  agreement  between  classical  and  quantum 
numerical  results  and  the  analytical  solution  given  by  eg.  (5a)  (  ). 

Fig.  11b  Probability  distribution  over  the  unperturbed  states  at  t  =  120 
for  the  case  of  fig.  I  la,  after  reversal  of  velocitiesa  at  t:=D0.  Notice 
that  the  quantum  system  (open  lozenges)  recovers  its  initial  state  to 
seventeen  digits  which  corresponds  to  numerical  errors.  In  contrast,  the 
classical  motion  (solid  lozenges)  proceeds  according  to  the  diffusion 
equation  (5a)  (squares). 

Fig.  12  Claasical  (solid  lozenges)  and  quantum  (open  lozenges)  ionization 
probability  (excitation  above  the  unperturbed  level  h:150)  as  a  function 
of  time  for  the  case  of  fig  1 1b.  Notice  the  perfect  specular  simmetry  of 
the  quantum  curve  about  the  time  of  reversal  t=dO. 

Fig.  13  Excitation  probability  at  time  t=  60  as  a  function  of  the  field 
intensity  for  different  values  of  n^,  and  Wo.  W,(t"ci)  =  Wj*i(eo)/W,*^i(€jj('))  is 
the  quantum  excitation  probability  at  Cq  rescaled  to  the  corresponding 
classical  excitation  probability  computed  at  Co=«q'''*.  'S  field 

the  rescaled  field  intensity.  {&.}  no=30,  <jOo=3;  (.)  no=45,  >'j.'o=1;  (  )  no=45, 
Wo=3;  (.)no=6D,  ho=66,  a>o=2;  (x)  no=66,  (jOo=3;  (o)  no=IOO,  Wo=3. 

The  fact  that  all  points  coresponding  to  different  no  and  Wq  meet  at  the 
value  W  =1  for  Cq^l  is  a  numerical  verification  of  our  estimate  (15);  it 
also  verifies  that,  in  the  delocalized  regime,  the  quantum  excitation 
probability  is  close  to  the  classical  value. 

Fig.  14a  Quantum  probability  distribution  TTn)  averaged  over  60  periods 
of  c  vrithin  the  interval  ci0'c:<l20.  Here  np^oD,  u'[|=2.  and  tp-O.O:- 
(ng.l4al);  t,;,=u.Oi:i(fig. I4a2);  eo-0.l4  (fig.  I4a3). 


Fiq.  Mb  5jme  Fiq.  Hj  with  n(,=  IOO,  'x'd-?.  Jnd  (Fiq.  Hh  I ); 

liiix'j.  I  inq.  (fig.  MtiJ). 

Fiq.  15  Excitjtiufi  prubdbilitq  Wj  ^  as  a  functiun  of  froquencq  'x'l-i  at 

r:qo  for  fi::ed  ty  -  0.04  and  for  different  iiy.  i.)  ny^TO:  ‘lj  ny  =  45:  i::f 
11,5=66;  (*)  0,5=100.  The  solid  line  gives  the  rlaoiiral  erritatinn 
probabilifq.  Notire  that,  by  iricreasinq  o.',,.  tne  quantum  excitation 
probability  becomes  much  less  than  the  correspundinq  classical  one  due 
to  the  fact  that  the  delocalization  border  (eq.  15)  increases  vvitti  '.i.),;,. 

Fig.  16  c,'uantum  probability  distribution  ftni  averaged  over  4'i  periods 
of  r  in  the  interval  40  ■r  <80  for  fixed  c, •,='!!. ciq,  n,',-,=i;i.7,  at,,]  different 
ny.  Fiq.  16a  ny=IOO  (_  .  ny=66  (-  Fig.  16b  0^=45  i_  .  nu=  i- 

The  classical  probability  distribution  is  also  shovvm--  •  In  order 

to  compare  the  quantum  distributions  with  different  n,-,.  with  the 
classical  one.  we  have  introduced  rescaled  quantities  f'  =  i  (1,5/66)  f  and  n 
=  i66/n(j)  n.  The  scaling  property  of  the  quantum  di stritiut ion  and  the 
fairly  good  agreement  with  the  classical  motion  is  due  t,:,  the 
delocal i zation  phenomenon. 

Fig.  17  Fine  structure  of  the  dependence  of  excitation  prpbabihty  '.v,  c; 
on  frequency  0.10  at  ~=40  for  fixed  cy=0.04.  io)  n,5=45:  (♦)  n,5=li>0:  (.) 
classical  results. 

Fig.  16  Ionization  probability  W  =2n-n  versus  field  frequency  (x),5 
after  a  time  r  -  qif.VY,  which  corresponds  to  the  same  real  physical  time 
t  for  all  frequencies.  We  have  set  nii^rjD,  n  =  'H9,  floreover, 

quantum  theory  classical  theory  io).  Nnti;e  that  'a',  is  here 

:.iimewhaf  less  than  n,|/z  becau;,e,  in  our  deiimtinn  of  thn  irmiration 
probability,  the  contribution  oi  states  with  n>n  is  also  included. 

Fig.  16  Same  as  fig.  16  'with  iiy=T.M.  ^|;,:■a.Cl75.  n  ='-»'. 
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Fig.  Llaii.icol  (.)  jurj  guafitutn  (♦,'  excitation  proDatnlitg  a;,  a  function 
of  time  for  n(,=66.  to^O.Oo.  olIq-  0.^3.  The  quantum  system  is  delocalized 
but  since  we  are  only  slightly  above  the  border,  the  quantum  excitation 
IS  less  than  the  classical  one. 


6.  Conclusions  and  recommendations 


riip  study  or  rue  i-dim  H  atom  in  a  nionochromatic  tield  that  we  have 
described  ih  the  present  paper  brings  into  the  light  a  hiirnber  of  facts  - 
some  of  w'hich  w'ere  rather  unexpected. 


These  facts  concern  both  the  actual  physics  of  atoms  in  microvvave 
fields,  and  the  yeneral  problem  of  quantum  dynamics  in  the  region  of 
classically  chaotic  motion.  Even  though  the  unperturbed  eigenfunctions, 
as  W'-eil  as  the  matrix  elements  oi  the  perturbation,  can  be  well 
approximated  by  tneir  semiclassical  expressions,  it  may  well  happen 
that  quantum  and  classical  time  evolutions  are  essentially  different,  due 
to  the  phenomenon  of  quantum  localization  of  chaos.  It  is  interesting  to 
note  that  investigations  of  this  phenomenon  were  prompted  by  studies  on 
the  rotator  model  /4/.  It  is  a  remarkable  fact  that  this  phenomenon, 
originally  detected  in  a  somewhat  artificial  model,  has  now  been  shown 
to  exist  in  a  physical  system,  so  that  there  is  a  real  possibility  to 
observe  it  in  laboratory  experiments  . 

On  the  other  hand,  for  the  H  atom  also  a  delocalization  regime  exists, 
and  our  theory  allows  determination  of  the  threshold  for  this  regime. 
Above  this  threshold,  the  excitation  of  the  guantum  system  can  be 
approximately  described  by  the  classical  diffusive  excitation.  This 
regime  of  excitation  is  much  more  efficient  than  the  direct  1 -photon 
ionization:  therefore  a  new  frequency  threshold  for  the  photoelectric 
eifect  appears,  which  is  determined  by  classical  border  for  frequency, 
■■0,;,  >  Actually  there  are  two  different  frequency  thresholds  and  co, 
.  so  that  strong  ionization  occurs  only  for  co^  The  latter 

threshold  ''.0|  IS  due  to  quantum  localization  of  classical  chaos. 


rne  delocalization  pnenomenon  explains  the  partial  success  oi 
classical  computations  in  reproducing  experimental  results  on  rnicrowa'-'e 
lOMization.  At  the  same  time,  how'ever.  the  localization  phenomerion  sets 
definite  limits  to  the  applicability  of  classical  mcdels,  which  are  due  to 


u  •  j  n  t  u  I M  lo  C  a  1 1  r  a  1 1  oil . 
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Although  3  discussion  oi  the  two-dim.  case  was  given  in  sec.  2.2, 
the  bulk  ot  the  results  presented  in  this  paper  were  related  to  the 
one-dimensional  case.  While  this  fact  does  not  certainly  affect  their 
conceptual  importance,  it  enforces  some  caution  when  comparing  them 
with  experiments  hiterto  performed.  Indeed,  an  analysis  of  the 
experiments  described  in  /4?/  shows  that  a  three-dimensional  theory  is 
required  to  model  them  properly.  A  different  series  of  experimerits/23/ 
is  amenable  to  a  one-dimensional  description,  but  here  a  static  electric 
field  CjCollinear  with  the  microwave  field  was  present  during  the  whole 

time  of  the  experiment.  In  order  to  correctly  model  these  experiments 
we  ought  to  add  a  term  zc^  in  our  Hamiltonian.  The  numerical  study  of 

the  resulting  quantum  dynamics  presents  some  technical  difficulty, 
because  it  seems  to  call  for  a  finer  description  of  continuous  spectrum 
than  allowed  by  the  method  exploited  here.  Perhaps  recourse  to  more 
sophisticated  numerical  methods  will  prove  necessary  but,  in  our 
opinion,  the  basic  qualitative  picture  of  localization-delocalization  will 
not  change. 

[lelocalization  is  also  a  challenging  subject  for  future  theoretical 
analysis.  This  phenomenon  has  been  predicted  on  the  grounds  of 
semiclassical  arguments,  which  are  best  suited  to  make  contact  with 
classical  chaotic  behaviour.  Nevertheless,  it  should  be  possible  to 
understand  it  in  purely  quantum  terms.  A  first  step  in  this  direction  may 
be  provided  by  a  recent  result/49/  that  a  qualitative  change  occurs  in 
the  numerically  computed  quasi-energy  eigenfunctions  of  the  one-dim. 
problem. 

A  few  concluding  remarks  are  in  order  concerning  the  relationship 
of  results  described  in  this  paper  to  the  general  themes  of  quantum 
chaos.  As  we  have  seen,  diffusive  excitation  and  ionization  are  brought 
about  in  the  classical  hydrogen  atom  by  the  onset  of  dynamical  chaos, 
which  is  a  regime  of  extreme  instability  of  traiectories  of  the  electron. 
A  physically  relevant  question  that  we  have  answered  above,  is  whether 
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and  so  on  -  survive  also  in  the  quantum  domain.  However,  the  more 
speculative  question  may  be  posed,  whether  also  arnjthing  oi  the 
conceptual  setup  ot  classical  chaos  -  instability.  irreversiDility,  and  so 
on  -  can  be  translated  in  a  quantum  context.  An  illustration  was  given  in 
this  paper  (sec.  3.5)  that  this  is  not  the  case.  However  similar  the 
quantum  evolution  may  appear  to  the  classical  (  insofar  as  the  population 
of  levels  is  concerned)  it  remains  strongly  stable,  in  sharp  contrast  to 
the  latter.  Therefore,  even  though  classical  chaos  was  shown  to  be 
relevant  in  predictin.,  lie  response  of  a  quantum  hydrogen  atom  to  an 
external  microwave  field,  it  must  be  stressed  again  that,  strictly 
speaKing,  no  true  chaos  is  possible  in  quantum  mechanics. 

We  shall  now  discuss  how  the  results  of  our  investigations  modify 
the  general  picture  of  quantum  chaos,  as  we  have  sketched  it  in  sec.  I. 
Indeed,  we  believe  that  in  the  light  of  these  findings  of  ours  some 
previously  accepted  views  must  now  be  modified;  at  the  same  time,  new 
developments  appear  now  possible,  which  hopefully  will  shed  light  on 
some  as  yet  scarcely  understood  phenomena  in  atomic  physics. 

Our  study  about  the  possibility  of  chaotic  diffusion  in  quantum 
periodically  perturbed  systems  led  us  to  definite  predictions  on  the 
existence  of  a  localization-delocalization  mechanism.  Further  theoretical 
analysis  and  experimental  work  is  required  in  order  that  this  mechanism 
may  enter  the  domain  of  physically  ascertained  facts,  and  indeed  work  is 
in  progress  in  both  directions.  I'leanwhile,  we  wish  to  stress  that  our 
theoretical  views  yield  considerable  clarification  in  an  otherwise  very 
confused  state  of  affairs.  Indeed,  judging  by  current  literature  and 
recent  international  meetings,  scientists  working  on  time-dependent 
problems  can  be  roughly  devided  in  two  categories.  In  tne  first  of  these 
we  classify  those  with  a  more  or  less  negative  attitude  in  regard  of  the 
possibility  that  some  cluotii,.  effects  may  survive  in  Quantum  Mecliamcs. 
It  includes  mainly  people  that  at  some  stage  of  their  scientific  work 
became  acquainted  with  the  kicked  rotator,  or  vvith  related  models,  and 
I'lave  therefore  a  direct  experience  of  the  stubborn  resistence  offered  by 
this  quantum  object  to  any  attempt  at  introducing  chaos  in  it.  These 
people  usually  develop  a  theoretical  attitude  according  to  which  the 


i]uantuiTi  suppression  of  dLinornicdl  chaos  is  not  lust  an  artefact  of  a 
queer,  highly  non-generic  model  but  has  deep  roots  in  the  very 
foundations  of  quantum  I'lechamcs.  Typical  representatives  of  this  class 
were  Hogg  and  Huberman,  who  even  telt  justified  in  dilating  the  rotor's 
stability  into  a  general  law,  according  to  which  any  quantum  system 
subject  to  a  periodic  perturbation  would  exhibit  a  strongly  recurrent 
behaviour,  with  the  only  possible  exception  of  non-generic  resonant 
situations. 

Even  though  this  contention  of  Hogg  and  Huberman  (71  turned  out  to  be 
an  erroneus  one  -  we  were  indeed  able  to  prove  that  the  rotor  itself  can 
be  nonrecurrent  without  being  resonant  (101  -  nevertheless  it  prompted 
important  developments.  The  claimed  impossibility  that  periodically 
perturbed  quantum  systems  may  follow  classical  chaotic  patterns  was 
assimilated  by  Fishman,  Grempel  and  Prange  [8]  to  another  more  firmly 
established  quantum  impossibility,  namely,  that  a  quantum  particle  in  a 
disordered  static  potential  may  ever  escape  to  infinity.  As  a  matter  of 
fact,  it  is  well  known  that  such  a  particle  will  stay  localized,  due  to  a 
complicated  interference  effect.  Fishman,  Grempel  and  Prange  were  able 
to  establish  a  formal  connection  between  time  dependent  problerris  and 
localization  problems  on  disordered  lattices,  that  proved  very  useful. 

It  is  Interesting  to  note  that  people  in  this  first  class  who  entered 
the  H-atom  problem,  at  least  initially  were  rather  skeptical  about  the 
possibility  that  anything  like  a  chaotic  ionization  may  exist  in  the  real 
quantum  H  atom,  and  considered  It  almost  obvious  that  some  more  or 
less  severe  limitations  would  be  imposed  by  quantum  effects  on  the 
chaotic  diffusion.  (Blumel  and  Smilanski,  (361  5hepe1yansky.  (11). 

A  somewhat  specular  attitude  is  displayed  by  workers  who  enter 
Quantum  Chaos  just  because  of  their  involvment  in  microwave  ionization. 
Being  aware  of  the  partial  success  of  classical  computatiuns  in 
reproducing  experimental  data,  and  relying  on  the  correspondence 
principle,  they  hardly  realize  that  a  theoretical  scheme  originally 
introduced  for  the  exotic  kicked  rotor  provides  elements  also  for  a 
quantum  theory  of  the  H  atom. 

According  to  our  results,  both  positions  ai.?  partially  justified.  The 
localization- delocalization  mechanism  provides  indeed  a  fey  to 
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leave  room  for  seemingly  diffusive  exitation  processes  even  in  quantum 
dynamics.  Much  work  is  now  required  in  order  to  give  firm  ground  to  our 
theoretical  results.  In  the  first  place,  our  theory  of  the  delocalization 
phenomenon  is  essentially  a  serniclassical  one.  It  would  be  nice  to 
understand  this  phenomenon  m  purely  quantum  terms,  i.e.,  by  making  as 
little  use  as  possible  of  the  underlying  classical  diffusive  picture.  In 
this  way  the  possibility  could  be  investigated,  that  a  similar  quantum 
instability  occurs  in  different  situations  for  which  no  classical  analog  is 
available-e.q..  in  problems  of  nuclear  physics. 

The  best  approach  towards  a  purely  quantum  theory  of  delocalization 
IS  ihe  spectral  analysis  of  the  quasi-energy  eigenvalues  and 
eigenfunctions.  Interesting  results  in  this  sense  were  recently  obtained 
by  Blumel  and  Srnilansky  (481  who  performed  a  numerical  computation  of 
quasi-energy  eigenfunctions  for  several  field  parameters  and  were  able 
to  detect  an  abrupt  change  in  their  shape  across  some  threshold.  We  are 
currently  investigating  the  possible  connnections  between  this 
phenomenon  and  delocalization. 

Another  important  question  is  hovs'  our  theory  should  be  modified  in 
order  to  apply  also  in  the  realistic  case  of  a  T>  dimensional  atom.  A 
3-dirnensions  numerically  solvable  quantum  model  is  then  needed, 
allowing  for  a  careful  analysis  of  the  localization  phenomenon  in  3 
dirnensionas.  We  have  already  some  theoretical  estimates,  according  to 
which  the  qualitative  picture  should  not  be  different  from  the  1  dim. 
case.  W'e  are  also  working  on  a  numerical  model,  which  is  still,  however, 
111  a  preliminary  stage. 

Insofar  as  our  theory  justifies  the  use  of  classical  concepts  under 
appropriate  conditions,  it  can  be  considered  as  a  theoretical  explanation 
of  the  previously  detected  phenomenon  ol  underthreshold  ionization,  lui 
the  other  hand,  the  loc alization  ptienornenon  in  the  hydrogen  atom  in  a 
microwave  field  is  a  re  lev. ant  prediction  of  ours  that  has  not  yet  been 
subjected  to  experimental  tests.  Moreover,  according  to  our  views,  .a 
scan  of  the  |■lnlzatlon  as  a  function  of  microwave  freguenoy  should 
e::po;.e  .3  big,  sharply  defined  peak  in  the  region  of  low  frequency  i  this, 
also,  calls  I  or  exfierimeni.al  veniic  ation. 


In  order  thot  a  detoiled  loborotord  checl.  Cjii  be  tnode,  we  i.liould 
include  in  our  model  some  minor  modifications  accounting  for  certain 
details  of  the  actual  experimental  setup,  for  example,  in  actual 
experiments  hydrogen  atoms  are  suDmitteo  to  the  combined  action  ot  the 
microwave  field  and  of  a  static  electric  field.  Even  though  the  latter 
should  not  essentially  modify  the  localization-delocalization  picture, 
Including  it  into  the  numerical  model  requires  great  care  and  calls  for 
extreme  computer  performances. 

In  drawing  a  final  balance  of  a  three-years  activity,  we  do  not 
certainly  feel  like  stating  that  the  major  problems  m  the  solution  ot 
which  we  purported  to  contribute  have  been  given  a  final  ansv/er 
Whatever. 

Nevertheless,  we  believe  that  the  world-wide  research  activity  on  the 
role  of  chaos  in  Quantum  Mechanics,  in  'which  'we  have  been  contributing, 
is  now  entering  an  almost  unexplored  realm  of  microphysics,  of  great 
potential  relevance  and  immediate  physical  interest. 
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8.  Appendixes 

Appendix  1  :  Quasi  Classical  riatnx  Elements  in  Paradcihc:  Coordinates. 


We  shall  here  get  the  expression  for  the  ^  coordinate  in  parabolic 
action-angle  variables  (n,,  112,  m,  X],  X^,  "I*).  To  this  end  we  introduce 
the  parabolic  coordinates 

x-  cosf 

siri'l'  (I.l) 

In  these  coordinates  the  unperturbed  Hamiltonian  takes  the  form; 

H=  2C/(^*T|)  p^^  ♦2'n/(^*Ti)  p.^.|‘'  *  I/(2^ti)  p^/  -2/(i;*T|.)  (T.2) 

The  transformation  to  action-angle  variables  (rii.n^.rrijXiA^.'j')  Is 
achieved  by  separation  of  variables  in  the  Hamilton  -Jacobi  equation  for 
■which  'we  refer  to  standard  textbooks  (  see  e.g.  30).  Here  v/e  just  recall 
that  the  generating  function  of  this  transformation  is  found  to  be: 

^  T| 

S(n„n2,rn,5,T(,'f)=  j  d^'  *  dri'  -rnT'  (1.3) 

where  the  canonical  momenta  Pj,  ,  P^are  given  by/30/: 


Pr  =  [E/2+  f;,/C  -  rrP/4i-] 

P^i  =  lE/2»  f.j/q  -m-/4iq-  1'''- 
P|;  ^  m 

P  =  '•’'1,2''  I  I 


E  :  -  l/  2rr 


n=  n,  ♦  n^  *  j  m 
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Then  the  angle  variables  X,,  Xj 

X|  2  “  ^5/^ni^2 


(1.4) 


can  be  obtained  by  differentiating  (1.3)  and  computing  the  integrals.  The 
procedure  is  greatly  simplified  by  the  introduction  of  the  auxiliary 
angles  Xi.X:  defined  by 


i  -  -2n^  p,  sin  Xi  +  2n  (n,  ♦  |m|/2) 
T|  =  -2n^'  P2  S'fi  Xa  2n  (02  »  |  m  |  /2). 


(1.5) 


where  the  parameters  Ui.U2  are  given  by: 
Mi.j  =  [h|,2  '>ri  -  n2,,)/n‘']'''^ 


In  this  way  we  get  the  following  result: 

X,  :  -  U|  cos  Xi  -  U2  cos  X2  ■  Xi 
X2  =  -  Pt  cos  Xi  -  ^2  cos  X2  -  X: 


fl.6) 


From  (1.5)  we  get  the  following  expression  for  z: 

’  =  1/2  iX’-g)  =  n-’tp^  sin  X2  -  Ui  s'n  Xi^  *  0  (n,  -  no)  (I.?) 


The  Luordinate  z  tan  be  expanded  in  a  double  Fourier  series  in  the 
angles  XiAz  coefficients  given  by: 

A,t,,  =  |d>-i  I'dx,  z  e 
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'.Ve  MOW  iutiititute  il.7)  for  2  in  this  intHijral,  jinj  chango  integration 
variables  to  Xi>  X;  tiy  using  (1.6).  Thus  we  find: 

2ti  in 

1,^  =  n‘  [  dx,  [  dXz  C'^Xi  fPz  sin  Xz  ‘  di  sin  Xt)  e  ''2^2'  ♦ 

0  0 

*  n  (n,  -  n,)  (].6) 

where  D(X,,  x^X)  =  d(X,.X;,.)/d(Xi.X2)  =  '  '  Mi  sin  Xi  ‘  M;.>  sin  x?  is  the 
Jacobian  determinant  for  the  transformation  (Xi,  x?)  (Xti  Xi'X 

Evaluating  the  double  integral  in  (1.8)  yields  formulas  (Id)  in  the  text. 

Appendix  11:  Solution  of  the  Fokker-Planck  equation. 

(n  order  to  solve  Eq.  (4)  'with  the  boundary  condition  <tf/ciii|^.-  =0, 

we  shall  first  perform  some  change  of  variables.  First  of  all,  putting 
t  2  .  yMi/'ni,.  the  Fokker-Planck  equation  tal  es  the  form: 

jf(y,t’)/dt  :  d/dij  (y^df/rJy) 

Now  let's  change  again  variables  to  z  ^  y''^'^  and  let's  introduce  a  new 
function  gfz,^)  according  to  f-z^  g.  This  (unction  q  must  then  satisfy: 

jq/jT  2  (1/4)  •>-q/dz-  +  ( l/4z)dg/7iz  -  ( l/z-')  g 

and  its  Laplace  rransiorm  g  '.r.s)  rnust  •satisiy  the  eguatiMn; 

.  (i,';)  /iq/.i-  -  (-4/-'  •4si  g  -  4g'.,sr'. 

A  Iiirfhei  1  timiji'  n!  'ar  i.irdf.  In  ■  =  .  ijudg,: 
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■j-q. ♦  (l/x)  jq/dx  -  ( I ’-t/x-','  q  -  -  (l/s)  q(x,0)  (II.  I ) 

By  the  same  changes  of  variables  we  find  that,  in  order  that  f  satisfies 
the  boundary  condition  df/dn|^.-  -0,  g  must  satisfy: 

3q(x,s)/3x ( -  -(2/x)  g(x^  s)  ,  where  2(s  n^/n)''^  (II. 2) 

The  general  integral  of  eq.  (II.  1)  can  be  written  as 

y  -  A  Iji'x)  ♦  B  K'?(x)  ♦  q 

Where  I2,  are . A  .B  are  numerical  constants  and  ^  is  a  particular 

integral  that  can  be  determined  e.  g.  by  Lagrange's  method: 

X 

y  :I;a'x)  [a  -  il/s)  J  x’g(x',0)  K2(x')dx'l  ♦ 

n 

X  (II. 3) 

♦  K,(X)  (B  ♦  (1/s)  [  X'g(X',0)  l2(X')ijX' 

0 


The  constants  A.B  can  then  be  chosen  so  that  the  boundary  condition 
(11.2)  is  satisfied.  Indeed,  upon  substituting  (II. 3)  into  (II. 2)  we  get: 


g(x,3)  =[  K,(x)  Ijixl/Iiix)]  ( 1/s)  f g(x',0)l2(x')x'dx'  ♦ 

_  0 

I^.(xi  (l/si  I  g(x',u)l  ^(x  ix  dx'*  l.;.(xi  (l/s)  ( q''.;-;'.Oil>(x')x’dx' 


(II.-l) 


Since  tin, Cl)  -•..ui'n,,),  we  must  cnooc-e  gcx.ci,'  in  nie  lorni  OHs/X'  - 1 ), 
Then  the  asymptotics  of  'll. 4'  for  s  ->  ■>!  and  fixed  y  has  the  form: 

g'-:.s'  i  4  i '■ '*1  1  V  - 1 1 1' -pi  s  I*'-'  .’v  -  v'-! ']  * 


whence  it  lollows 


exp[2 /3(x/2 /s- 1  )1  ( 

f(y,s)  "  1/(2  ys)  (exp[2v's(|t|/Vy-2/-v''^l'* 

b->  00 

+exp[2y  s()/V'y- 1)1 1 

£q.  (5a)  in  the  text  easily  follows  from  the  last  formula. 

Appendix  III:  Estimate  for  the  Delocalization  Border. 

Let's  evaluate  the  2nd  rnornent  of  the  distribution  over  the  levels: 
i'll  =  <(An)“>  =  ^■(n-<n>)‘>.  From  the  diffusion  equation  w'e  get 
approximately 

d/dt  <(,Ah)‘'>  %  <D>=  arig^<(,n/nf|)^  >  (Ill.l) 

where  a=  2  The  equation  for  the  1st  moment  gives 

d<h>/dt  :  (3  3  no/2)  <(n/no)^>  (III. 2) 

In  order  to  solve  (lII.l)  and  (in.2)  we  will  use  a  rough  approximation, 
namely,  we  will  substitute  for  n  its  mean  value  <n>.  Doing  so,  and 
performing  the  integration,  we  obtain 

-ri>  -  riol  I- 3a t/2r' 

<t  An)‘>=  riiyif  1  -~ar  /2)’^- 1 1/3 

Tlio  localization  condition  ■  o'  -  (An)-'  ~  o  l'  gives  an  equation  lor 
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t,  the  least  coot  of  which  deterinines  the  local Ization  lenqth  1: 

tno‘/3)  !( 1  -3ta/2)'2  -11  =  tV(x‘  (III.3) 


1  =  r/(x 


A  straigthforward  manipulation  gives  then  formulas  (13),  (14). 


Appendix  IV:  A  Method  for  computing  Hypergeometric  functions. 

The  numerical  computations  of  matrix  elements  in  (32)  ?  presents 
some  technical  difficulty  since  a  direct  expansion  of  hypergeometric 
functions  in  series  of  powers  of  -4nno/(n-n(i)^  doesn't  give  correct 
values  of  for  n~s  ~  100  due  strong  cancellations  of  different  terms 

and  finite  computer  precision.  Therefore  in  order  to  compute  the 
hypergeometric  function  F  we  used  a  different  method  based  on  the 
recursion  formulas  between  values  of  F  for  three  consecutive  values  s-l, 
s,  s*-l  (  see,  e.g./39/).  The  method  is  essentially  as  follows: 

We  take  two  values  of  F  for  s=0,  s=l  and  then  we  recurrently 
determine  all  F  up  to  $=0^  ss  n  .  After  that  we  take  two  arbitrary  values 
for  F^^  and  F^^,  for  s  =m  and  s=m+l  where  m  was  chosen  (for 

example  m  ss  5  s^^^  )  and  recurrently  determined  values  of  F  for  nj-iO< 
s  rn  .  These  latter  values  of  F  differ  from  the  actual  ones  only  by  a 
numerical  factor  Cf.  The  value  of  this  constant  was  obtained  by 
comparison  with  F  computed  for  s<  n,.  For  different  values  s<  n.^  the 
constant  Cp  was  obtained  with  a  precision  s;  )0’’“  After  taking  into 
accouiit  this  coiistaiil  factor  we  got  precise  values  of  F  for  '7 
which  didn’t  change  upon  changing  the  arbitrary  values  F^,  .  F^^,,  . 
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